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The  research  described  herein  was  conducted  by  Mechanical 
Technology  Incorporated  under  the  terms  of  Contract 
DAAJ02-72-C-0040.  The  work  was  performed  under  the  technical 
management  of  Mr.  Leroy  T.  Burrows,  Technology  Applications 
Division,  Eustis  Dir"  'torate. 

During  the  past  decade,  vibration  and  noise  measurements 
and  data  reduction  procedures  have  improved  to  the  point 
where  it  can  be  (and  has  been)  clearly  shown  that  noise 
and  vibration  are  directly  relatable  to  each  other.  Moreover, 
many  noise  components  have  been  shown  to  be  directly  relatable 
to  the  gear  mesh  frequencies  in  such  drive  trains,  and 
analytical  methods  have  been  formulated  to  help  understand 
and  control  them.  These  methods  deal  with  the  mechanical 
vibrations  of  the  gearbox  components.  Other  significant 
signals  that  are  present,  however,  are  not  directly  relatable 
to  the  mesh  frequencies.  Some  of  these  signals,  called 
"sidebands",  have  been  found  to  occur  in  tests  of  helicopter 
rotor-drive  gearboxes. 

The  major  aims  of  this  study  were  as  follows: 

1.  Using  existing  measured  and  calculated  CH-47  ring 
gear  acceleration  data  together  with  sideband 
amplitude  prediction  methods  presently  under 
development,  the  contractor  was  (a)  to  investigate 
CH-47  lower  planetary  mesh  planet-pasB  sideband 
amplitudes  for  several  lower  plane t-to-ring 

gear  mesh  relationships  in  order  to  determine 
sensitivity  of  sidebands  to  planet  phasing,  (b)  to 
identify  other  design  parameters  expected  to  be 
useful  in  controlling  planet-pass  sidebands  and 
describe  their  effects  and  importance,  and  (c)  to 
dray  conclusions  concerning  mechanisms  producing 
planet-pass  sidebands. 

2.  The  computer  program  entitled  GEARO,  which  is  in 
the  possession  of  the  contractor  and  the  Government, 
was  to  be  modified  »nd  extended. 

Appropriate  technical  personnel  of  this  directorate  have 
reviewed  this  report  and  concur  wi'Ji  the  conclusions  contained 
herein. 
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SUMMARY 


Users  of  geared  power  trains  have  begun  to  recognize  the  importance  of  the 
high-frequency  vibrations  which  are  present  in  virtually  all  operating  gear¬ 
boxes:  these  vibrations  are  the  key  to  understanding  the  gearbox  condition 
in  real  time,  a  problem  of  considerable  current  importance.  An  immediate 
outward  indication  of  the  presence  of  such  vibrations  is  the  noise  produced 
by  a  gearbox.  Even  the  untrained  ear  can  distinguish  the  presence  of 
signals  which  arise  in  the  gearbox.  (Sensors  are  of  coutse  required  to 
obtain  information  of  the  quality  required  for  engineering  purposes.) 

During  the  past  decade,  vibration  and  noise  measurements  and  data  reduction 
procedures  have  improved  to  the  point  where  it  can  be  (and  has  been)  clearly 
shown  that  noise  and  vibration  are  directly  relatable  to  each  other.  More¬ 
over,  many  noise  components  have  been  shown  to  be  directly  relatable  to  the 
gear  mesh  frequencies  in  such  drive  trains,  and  analytical  methods  have  been 
formulated  to  help  understand  and  control  them.  These  methods  deal  with  the 
mechanical  vibrations  of  the  gearbox  components. 

Other  significant  signals  which  are  present,  however,  are  not  directly  re¬ 
latable  to  the  mesh  frequencies.  Some  of  these  signals,  called  "sidebands", 
have  been  found  to  occur  in  tests  of  helicopter  rotor-drive  gearboxes.  As 
originally  conceived,  this  investigation  had  the  rather  limited  objective 
of  providing  the  designer  of  geared  power  trains  with  an  analytical  tool 
which  could  be  used  to  predict,  and  thus  control,  the  frequencies  and 
amplitudes  at  which  vibration  sidebands  are  produced  by  operating  gearboxes. 
As  the  work  progressed,  however,  it  became  apparent  that  the  analyses  and 
associated  understanding  could  also  have  a  significant  and  far-reaching 
impact  on  the  more  general  problem  of  on-line  monitoring.  It  is  important, 
therefore,  that  the  results  of  this  study  be  viewed  in  the  context  of  their 
potential  impact  on  this  technology  area,  as  well  as  upon  the  area  of  gear¬ 
box  noise  reduction. 

The  major  aims  of  this  study  were  as  follows: 

1.  Development  of  an  engineering  understanding  of,  and  methods  for 
predicting,  geometrically-induced  planet-pass  vibration  sidebands 
which  accompany  normal  planetary  gear  reduction  operation;  and 

2.  Development  of  an  engineering  understanding  of,  and  new  analytical 
methods  for  treating,  the  vibration  sidebands  which  are  produced 
by  undesirable  gear  characteristics,  such  as  tooth  support  dis¬ 
continuities  (cracks),  gear  runout,  and  variations  in  tooth 
transmitted  forces . 

These  objectives  have  been  achieved.  Sidebands  as  a  second  major  category 
of  geared  drive  train  vibration  signals  can  now  be  described  and  discussed 
directly  m  terms  of  hardware  condition.  An  engineering  understanding  of 
the  mechanisms  which  cause  many  of  the  kinds  of  vibrations  produced  by 
gear  meshes  thus  exists . 

The  payoff 8  from  such  an  understanding  can  be  enormous.  First,  gear  train 


designers  now  have  important  vibration  analysis  tools  for  minimizing,  at  tne 
time  of  design,  the  disturbances  produced  by  gear  meshes.  This  will  not 
only  make  gear  trains  quieter,  but  till  also  reduce  their  internal  forces, 
with  significant  improvements  in  lifetimes  of  all  components. 

Second,  when  properly  exploited  by  joint  technologist/manufacturer  teams, 
the  high-frequency  vibration  analysis  tools  will  permit  real-time  condition 
monitoring  to  be  approached  on  an  individual  signal  component  basis  in  which 
the  precise  meaning  of  each  component  is  well  understood,  rather  than  on 
the  multiple  component  basis  associated  with  usual  signature  analysis 
techniques.  The  engineering  application  of  the  vibration  analysis  tools 
permits  specific  signal  components  to  be  explained  on  a  detailed  basis. 

This  removes  the  uncertainty  associated  with  the  use  of  signal  level  ratio 
techniques  and  the  need  for  lengthy  test-bed  study  programs  involving 
failure  implants. 

Extension  of  the  gear  mesh  analysis  techniques  to  high-contact-ratio  gearing 
is  now  in  order,  as  is  an  extension  of  the  torsional  response  analysis  to 
coupled  torsional-lateral-axial  vibrations. 
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Dr.  Robert  H.  Badgley  of  Mechanical  Technology  Incorporated  served  as  Pro¬ 
gram  Manager  for  the  efforts  reported  herein.  The  contract  was  carried  out 
under  the  technical  cognizance  of  Mr.  R.  Burrows,  Eustis  Directorate,  U.  S. 
Army  Air  Mobility  Research  and  Development  Laboratory,  Fort  Eustis,  Virginia. 

Special  credit  is  due  to  Mrs.  L.  Cziglenyi  of  MTI,  who  programmed  the  com¬ 
puter  program  modifications  and  who  conducted  the  calculations. 


TABLE  OF  CONTENTS 


SUMMARY .  iii 

FOREWORD  .  v 

LIST  OF  ILLUSTRATIONS .  ix 

LIST  OF  TABLES .  xi 

INTRODUCTION  .  1 

DESCRIPTION  OF  PROGRAM .  3 

Task  I  -  Investigation  of  Planetary  Mesh  Planet-Pass 

Sidebands  .  3 

Task  II  -  Analysis  of  Gear -Mesh -Induced  High-Frequency 

Vibration  Spectra  and  Calculations  .  3 

INVESTIGATION  OF  PLANETARY  MESH  PLANET-PASS  SIDEBANDS .  4 

Analysis  of  Planet-Pass  Induced  Vibration  .  4 

Planet-Pass  Induced  Vibration  and  Sideband  Amplitude 

Calculations  .  9 

ANALYSIS  OF  GEAR-MESH -INDUCED  HIGH-FREQUENCY  VIBRATION  SPECTRA 

AND  CALCULATIONS .  20 

Analysis  of  Gear  Mesh  Excitation .  20 

CH-47  Spiral  Bevel  Gear  Mesh  Calculations  .  22 

CH-47  Lower  Planetary  Plane t-to-Ring  Gear  Mesh 

Calculations  . .  29 

CH-47  Spiral  Bevel  Gear  System  Response  Calculations  ....  34 

DISCUSSION  OF  RESULTS .  39 

Discussion  of  Task  I  Results .  39 

Discussion  of  Task  II  Results .  40 

CONCLUSIONS . 44 

RECOMMENDATIONS .  46 

LITERATURE  CITED  .  47 


vii 


APPENDIXES 


I  CALCULATION  OF  POWER  SPECTRAL  DENSITY  FUNCTION .  49 

II  COMPUTER  PROGRAM  FOR  PREDICTION  OF  GEAR  MESH  EXCITATION 

SPECTRA .  52 


DISTRIBUTION 


79 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1  CH-47  Forward  Rotor-Drive  Gearbox  Lower  Planetary 

Gear  System . . .  5 

2  Gear  System  In  a  Planet-Fixed  Frame  and  Corresponding 

Force  Diagram .  6 

3  Dynamic  Model  of  CH-47  Lower  Planetary  Ring-Gear  Casing  .  H 

4  Rinf-Gear  Casing  Response  Functions  .  12 

5  Radial  Vibration  Amplitude  vs.  Time  at  Fixed  Location 
on  CH-47  Forward  Rotor-Drive  Gearbox  Ring  Gear  for 

t  =  37.6° .  13 

6  Frequency  Spectra  for  i(i  »  0°  and  37.6° .  1^ 

7  Radial  Vibration  Amplitude  vs.  Time  at  Fixed  Location 
on  CH-47  Forward  Rotor-Drive  Gearbox  Ring  Gear  for 

I  =  0° .  16 

8  Radial  Vibration  Amplitude  vs.  Time  at  Fixed  Location 
on  CH-47  Forward  Rotor-Drive  Gearbox  Ring  Gear  for 

*  -  90° .  17 

9  Radial  Vibration  Amplitude  vs.  Time  at  Fixed  Location 
on  CH-47  Forward  Rotor-Drive  Gearbox  Ring  Gear  for 

f  -  142.4° .  18 

10  Frequency  Spectra  for  f  =  90°  and  142.4° .  1® 

11  Computer  Program  Structure .  23 

12  Frequency  vs.  Excitation  Amplitude  for  Pinion  Shaft 

Runout,  CH-47  Spiral  Bevel  Gear  Mesh .  26 

13  Frequency  vs.  Excitation  Amplitude  for  Pinion  Shaft 

Runout,  CH-47  Spiral  Bevel  Gear  Meth.  . .  27 

14  Frequency  vs.  Excitation  Amplitude  for  Tooth  Force 

Variation,  CH-47  Spiral  Bevel  Gear  Mesh  . .  28 

15  Frequency  vs.  Excitation  Amplitude  for  Planet  Runout, 

CH-47  Lower  Planetary  Planet-to-Ring  Gear  Mesh .  30 

16  Frequency  vs.  Excitation  Amplitude  for  Tooth  Force 
Variation,  CH-47  Lower  Planetary  Planet-to-Ring  Gear 

Mesh .  32 


lx 


LIST  OF  ILLUSTRATIONS  (CONCLUDED 


Figure 

Page 

17 

Frequency  vs  Excitation  Amplitude  for  Tooth  Support 
Compliance  Variation,  CH-47  Lower  Planetary  Planet-to- 
Ring  Gear  Mesh . 

.  33 

18 

Calculated  Dynamic  Tooth  Force  at  CH-47  Spiral  Bevel 

Gear  Mesh  . 

.  37 

19 

Calculated  Dynamic  Tooth  Force  at  CH-47  Planet-Ring 

Gear  Mesh . . . 

.  38 

20 

Vibration  Spectrum  Measured  Near  the  Spiral  Bevel  Gears 
(From  Reference  8) . 

.  42 

>» 


LIST  OF  TABLES 


Table 

I 

II 

III 


Conversion  of  Spiral-Bevel  Gears  Into  Equivalent 
Spur  Gears  for  CH-47  Forward  Rotor  Transmission 
(From  Reference  2) . 

CH-47  Lower  Planetary  Planet  and  Ring  Gear  Parameters  . 

CH-47  Spiral  Bevel  Mesh  Excitation  Amplitudes  and 
Corresponding  Predicted  Peak  Dynamic  Tooth  Forces  at 
Indicated  Frequencies  . 


24 

29 


35 


xi 


4 

3 


I 


INTRODUCTION 


Geared  power  train  vibrations  can  occur  at  many  different  frequencies.  The 
most  common  of  these  vibrations  may  be  found  at  the  mesh  frequencies,  and 
their  integer  multiples,  of  particular  gear  meshes  in  the  train.  These  are 
known  to  be  caused  by  the  mesh  properties  of  the  gear  teeth.  Less  well 
understood  are  vibrations  which  occur  at  the  foregoing  frequencies  plus  and 
minus  Integer  multiples  of  other  frequencies.  Such  components  are  called 
sidebands.  While  methods  had  been  developed  for  predicting  the  levels  of 
the  vibration  components  at  the  mesh  frequency  and  its  integer  multiples  [l 
through  7],*  these  methods  were  not  capable  of  treating  the  sideband  com¬ 
ponents.  Unfortunately,  high  acoustic  noise  levels  can  be  produced  by  both 
types  of  components. 

In  addition,  there  is  increasing  recogniticn  of  the  fact  that  vibrations 
which  produce  noise  also  carry  information  about  the  dynamic  behavior  of  the 
drive  train  [8  through  11],  in  effect,  information  from  which  the  condition 
of  the  drive  train  components  may  be  inferred.  Such  vibrations  are  known  to 
be  present  in  virtually  all  geared  systems,  and  they  have  been  recorded  and 
monitored  for  diagnostic  purposes  for  many  years  by  many  people.  However,  a 
detailed  engineering  understanding  of  the  meanings  of  the  shapes  and  ampli¬ 
tudes  of  the  measured  spectra  has  not  proceeded  in  company  with  the  develop¬ 
ment  of  methods  for  sensing  and  displaying  these  spectra. 

Sidebands  are  produced  during  normal  operation  of  gearboxes  which  employ 
planetary  reductions.  Such  sidebands  normally  occur  at  the  planetary  mesh 
frequency  plus  and  minus  the  planet-pass  frequency.  (They  may  also  be  sim¬ 
ilarly  distributed  about  twice  mesh  frequency,  three  times  mesh  frequency, 
etc.)  It  must  be  stressed  that  their  presence  is  due  simply  to  the  kinemat¬ 
ics  of  the  planetary  reduction,  wherein  the  planets  physically  pass  any 
stationary  point.  The  presence  of  each  planet  in  turn  changes  the  vibration 
properties  of  the  gearbox  structure,  both  from  an  impedance  viewpoint  and 
also  more  importantly  from  the  amount  of  excitation  applied  to  the  structure 
by  the  moving  mesh.  This  change  is  periodic  at  planet-pass  frequency. 

Sidebands  are  also  produced  during  normal  operation  of  gearboxes  with  gear 
meshes  at  more  than  one  frequency.  These  sidebands  are  normally  found  at 
one  mesh  frequency  plus  and  minis  integer  multiples  of  the  other  mesh  fre¬ 
quencies  (and  of  course  at  other  frequency  combinations  as  mentioned  above). 
Such  sidebands  are  caused  primarily  by  the  dynamic  properties  of  the  drive 
train  components  (i.e.,  coupled  torsional-lateral-axial  vibrations),  which 
permit  dynamic  tooth  force  variations  to  occur  in  one  mesh  at  frequencies 
of  other  meshes. 

Other  sources  of  sidebands  do,  of  course,  exiBt,  but  these  are  for  the  most 
part  associated  with  the  presence  of  undesirable  component  behavior.  Per¬ 
haps  the  best-known  sideband  source  is  that  due  to  runout  of  a  gear  because 
of  machining  or  assembly  inaccuracies.  This  type  of  sideband,  which  the 


*Number8  in  brackets  refer  to  literature  cited  at  the  end  of  this  report. 


analysis  described  herein  can  predict,  is  typically  found  at  mesh  frequency 
plus  and  minus  shaft  rotation  frequency.  It  can  be  produced  by  runout  not 
only  of  pinion  or  gear  in  a  simple  mesh,  but  also  of  a  planet  gear  relative 
to  its  bearings  on  a  planet  carrier. 

Other  undesirable  effects  can  also  produce  sidebands.  For  instance,  varia¬ 
tion  of  tooth  support  stiffness  around  the  circumference  of  a  gear  can  alter 
the  mesh  properties  in  a  manner  which  repeats  at  gear  running  speed,  and 
which  can  have  many  forms  depending  on  the  circumferential  distribution  of 
the  stiffness  variation.  A  typical  cause  of  such  sidebands  would  be  a 
cracked  gear  web,  which  the  analysis  described  herein  can  treat. 
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DESCRIPTION  OF  PROGRAM 


In  this  study,  f>e  characteristics  of  vibration  sidebands  produced  by  gear 
meshes  in  both  single  gear  meshes  and  planetary  gear  reductions  were  in¬ 
vestigated.  The  study  was  conducted  in  two  tasks,  the  details  of  which  are 
described  below. 

During  the  course  of  earlier  test  efforts,  it  became  obvious  that  gearbox 
vibration  and  noise  signal  components  were  being  produced  at  frequencies 
which  corresponded  to  those  at  which  various  sideband  signals  were  expected. 
It  was  recognized  that  these  signals  would  have  to  be  explained  analytical¬ 
ly  before  they  could  be  dealt  with  properly.  Hence  the  derision  was  made 
to  treat  the  signals  at  their  sources,  i.e.,  the  gear  mesh  itself. 

Since  earlier  analytical  efforts  had  yielded  a  computer- implemented  analysis 
for  predicting  the  vibration  excitation  properties  of  normal  gear  meshes, 
the  decision  was  made  to  upgrade  this  computer  program  to  incorporate  the 
new  analyses.  The  upgraded  computer  program  was  modularized  so  as  to  per¬ 
mit  future  inclusion  of  other  gear  types  and  effects.  A  schematic  diagram 
of  the  computer  program's  capabilities  is  presented  and  discussed  later  in 
the  report  (as  Figure  11).  The  analyses  described  herein  have  been  included 
in  this  program. 

TASK  I  ~  INVESTIGATION  OF  PLANETARY  MESH  PLANET-PASS  SIDEBANDS 

As  a  result  of  earlier  studies,  predicted  dynamic  behavior  of  the  CH-47  for¬ 
ward  rotor  drive  gearbox  ring  gear  was  available.  This  predicted  data  was 
used  to  study  the  complex  vibrations  existing  at  a  preselected  point  on  the 
ring  gear  (corresponding  to  a  location  where  measured  data  had  been  taken). 

These  studies  produced  radial  vibration  levels  versus  time,  and  correspond¬ 
ing  vibration  spectra,  for  various  planet  phasing  relationships  in  the  low¬ 
er  stage  planetary  reduction.  Planet  phasing  is  under  the  control  of  the 
gear  designer,  and  thus  it  can  be  altered  to  modify  and  thus  reduce  mesh 
frequency  vibration  sidebands. 

TASK  II  -  ANALYSIS  OF  GEAR-MESH- INDUCED  HIGH-FREQUENCY  VIBRATION  SPECTRA 
AND  CALCULATIONS 


This  phase  of  the  study  treated  both  the  spiral  bevel  and  lower  stage 
planetary  planet- to-ring  gear  meshes.  Spiral  bevel  gear  shaft  runout  and 
externally- imposed  tooth  mesh  force  variations  were  studied,  and  gear  mesh 
excitation  spectra  produced  by  these  effects  were  predicted. 

In  the  care  of  the  planet-to-ring  mesh,  planet  runout  and  externally- imposed 
tooth  mesh  force  variations  were  considered,  and  excitation  spectra  caused 
by  these  effects  were  predicted.  In  addition,  the  mesh  excitation  spectrum 
resulting  from  the  condition  where  a  number  of  consecutive  ring  gear  teeth 
have  relatively  soft  support  stiffness  (such  as  could  be  caused  by  a  local 
crack)  was  predicted. 
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INVESTIGATION  OF  PLANETARY  MESH  P1ANET-PASS  SIDEBANDS 


ANALYSIS  07  PLANET-PASS  INDUCED  VIBRATION 

A  schematic  of  the  planetary  gear  system  to  be  analyzed  is  shown  in  Figure 

1.  The  ring  gear  is  stationary  and  the  sun  gear  rotates  with  speed  . 

s 

The  four  planets  rotate  about  their  own  centers  with  a  speed  of  a>  (relative 

o 

to  the  planet  carrier),  and  the  planet  carrier  rotates  with  speed  (planet 

orbiting  speed).  Depending  on  the  numbers  of  teeth  on  the  component  gears, 
the  four  planets  are,  in  general,  not  equally  spaced.  In  the  particular 
case  considered  herein,  6  is  the  angular  offset  of  the  pair  (C,D)  with  re¬ 
spect  to  the  pair  (A,B). 

In  this  analysis  the  ring  gear  is  considered  as  an  elastic  shell-type  struc¬ 
ture.  At  a  given  location,  the  ring  gear  experiences  a  dynamic  tooth  mesh 
force  each  time  a  planet  passes.  This  mesh  force  is  associated  with  the 
tooth  mesh  frequency  resulting  from  the  transfer  of  load  from  one  pair  of 
teeth  to  the  other.  In  this  calculation,  it  is  assumed  that  the  dynamic 
tooth  mesh  force  has  a  rectangular  pulse  form  whose  magnitude  varies  sinus¬ 
oidally  with  the  mesh  frequency  (see  below).  This  oscillating  force  pro¬ 
duces  an  oscillating  deformation.  It  is  this  planet-pass  induced  vibration 
of  the  ring  gear  that  is  treated  in  this  analysis. 

In  a  coordinate  frame  fixed  with  the  planrt  carrier,  both  sun  and  ring  gears 
rotate  as  shown  in  the  upper  diagram  of  Figure  2.  The  sun  gear  rotates  with 
a  speed  of  —  o^,  while  the  ring  gear  moves  in  the  opposite  direction  with 

the  planet  orbiting  speed  u^.  In  this  reference  frame,  the  forces  acting  on 

the  ring  gear  due  to  the  planet-ring  gear  meshes  occur  at  fixed  angular 
locations  of  6  ■  0,  tt/2  -  6  ,  tt,  3tt/(2)  —  6,  corresponding  to  the  locations 
of  planets  A,  D,  B  and  C  respectively. 

Since  planets  A  and  B  are  in  phase,  the  dynamic  force  per  unit  area  due  to 
gear  meshes  at  these  two  locations  can  be  approximately  expressed  as 

Fab  (6 ,t)  -  PAB  (0)  cos  ^t  (1) 


where  t  is  the  time,  PAB  the  normal  pressure  due  to  normal  tooth  mesh  forces, 
and  f^  thf  tooth  mesh  frequency  defined  by 

Hi  m  ®o  "p  "  “b  “r  <2> 


In  the  above,  N  and  N^  are  the  numbers  of  teeth  on  the  planet  and  ring 


gear  respectively.  It  may  be  assumed  that  P  consists  of  two  identical 
rectangular  pulses  occurring  at  the  locations  of  planets  A  and  B,  and  that 
the  width  of  each  pulse  is  one  tooth  spacing,  as  shown  in  the  lower  diagram 
of  Figure  2.  The  magnitude  of  the  rectangular  pulse  P  is  related  to  com¬ 
mon  gear  parameters  as  follows:  0 
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Figure  1.  CH-47  Forward  Rotor-Drive  Gearbox 
Lower  Planetary  Gear  System. 


5 


vwww 


p  * 

0 


F  tan  0 


*v 

(3) 


where  F  is  the  tangential  tooth  force,  0  the  pressure  angle,  r^  the  pitch 
radius  of  ring  g;ar,  and  1.  th?  tooth  face  width. 

As  a  result  of  the  angular  iffset  b,  there  exists  a  temporul  phase  differ¬ 
ence  between  the  dynamic  forces  at  planets  C  and  D  and  those  at  A  and  B. 

Since  rotation  of  the  sun  gear  over  one  tooth  spacing  corresponds  to  one 
full  cycle  in  gear  meshing,  it  can  be  shown  that  the  temporal  phase  lag  of 
the  tooth  meshes  of  planets  C  and  D  relative  to  those  of  A  and  B  is 

¥  -  6  •  N  (4) 

8 

where  N  is  the  number  of  teeth  of  the  sun  gear.  The  dvnamic  forces  at  planets 
s 

C  and  D  can  be  thus  expressed  as 

FCD  <6-  '>  *  PCD  <9>  COS  (fM  t-’,)  (5> 

where  P  is  similar  to  P  and  is  shown  in  the  lower  diagram  of  Figure  2. 

CD  AB 

There  exist  radial  displacement  response  functions  for  the  ring  gear  shell 

due  to  the  dynamic  forces  F.  and  F  .  These  response  functions  depend  on 

AB  CD 

the  elastic  characteristics  of  the  shell-type  ring  gear  structure.  Neglect¬ 
ing  ring  gear  inertia  (valid  for  shell-type  structures),  the  response 
function  due  to  the  meshing  of  planets  A  and  B  is  in  phase  with  its  forcing 
function  F^fi.  It  may  in  general  be  written  as 

bAB  cos  Si  1  (6) 

where  z  is  the  axial  coordinate,  designated  for  the  general  case  that  the 
ring  gear  shell  is  nonuniform  axially.  Similarly,  the  response  function 
due  to  the  meshing  of  planets  C  and  D  is 

bCD  ^6,Z^  C0S  t_Y*  <7) 

The  total  response  function  for  the  ring  gear  shell  is  the  sum  of  Equa¬ 
tions  (6)  and  (7),  i.e., 

w  (6,z,t)  -  bAg  cos  ^  t+bCJ)  cos  t-Y)  (8) 

where  w  represents  the  dynamic  radial  response  at  a  location  whose  coordinates 
are  (0,z)  on  the  ring  gear  casing. 
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After  some  manipulations,  Equation  (8)  becomes 

w  (0,2, t)  ■  A  (0,z)  cos  fM  t+B  (0,z)  sin  fM  t  (9) 

where 

A  (0,z)  “  bAB  +  bCD  cos  Y  (10) 

B  (6,0  -  bCD  sin  Y  (11) 

If  the  ring  gear  is  axially  symmetric,  both  b.  and  b  are  even  function3 

in  0  and  have  a  periodicity  of  n  in  0,  and  so  do  the  functions  A  (6,z)  and 
B  (@,2).  Therefore,  the  latter  may  be  expanded  in  terms  of  Fourier  series, 
i.e., 


and 


A  (0,z) 


6sl— 
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CO 
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E 

n*l 


2n 


(z)  cos  2n  0 


B  (0 ,z) 


bo  <z> 
2 


+  E 
n-1 


b2n  (z)  cos  2n  0 


(12) 


(13) 


Let  0q  be  the  angular  coordinate  of  a  fixed  point  on  the  ring  gear.  Since 
the  ring  gear  rotates  clockwise  with  speed  cu^  (see  Figure  2), 

0q  (t)  -  0t  +  i-  (14) 


where  0^  is  the  coordinate  of  the  fixed  point  at  t  *>  0.  Setting  0^  =  0, 

0q  (t)  -  a^t  (15) 


Thus,  the  fixed  point  on  the  ring  gear  travels  clockwise  starting  from  the 
position  of  8  ■  0  in  the  frame  fixed  with  the  planet  carrier. 

To  calculate  the  vibration  at  a  fixed  point  on  the  ring  gear,  Equation  (15) 
is  substituted  for  0  into  the  total  response  function,  Equation  (9).  This 
yields 


wq  (t)  -  w  (o^t,  zq,  t) 


A  (u^t. 


zo)  cos  ^  t 

+  B  (tt^t,  zq)  sin  fM 


(16) 

t 


where  zq  is  the  axial  coordinate  of  the  fixed  point. 

Using  the  Fourier  representations  of  Equations  (12)  and  (13), 
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2  a2n  [C°8  (fM  +  2n  V  1  +  COS  (fM  “  2n  V  *] 

2  JL  b2n  [8i0  (fM  +  2n  V  C  +  8in  (fM  "  2n  V  fc] 


(17) 


The  above  equation  represents  the  planet-pass  induced  vibration  amplitude 
history  at  a  fixed  point  on  the  ring  gear.  At  t  =  0,  planet  A  is  at  the 
location  of  this  fixed  point;  and  as  time  goes  on,  planets  D,  B,  and  C  pass 
the  point  in  sequence. 


From  Equation  (17),  the  vibration  amplitude  at  gear  tooth  mesh  frequency  is 

12  2  ^/2 

M  «  f  <•  +  bo  >  U8) 

O  L  0  O 

The  amplitudes  at  various  sidebands  about  the  mesh  frequency  are 

1  2  2  1/2 

Mn  "  2  (a2n  +  b2n  >  for  n  "  1.2,3...  (19) 

Since  there  are  four  planets,  the  planet-pass  frequency  is 

fp  "  4  %  (20) 

f 

The  sidebands  occur  at  +  n  with  n  ■  1,2,3,....  In  the  case  of  zero 

offset  (6  ■  0),  sidebands  would  appear  only  at  f^  +  n  f ^ ,  because 

the  function  A  and  B  would  have  a  periodicity  of  tt/2  instead  of  tt.  Physi¬ 
cally  it  means  that  at  the  fixed  observation  point,  the  pass  of  planets  C 
and  D  is  Identical  to  that  of  planets  A  and  B. 


Equation  (17)  gives  the  amplitude-time  relationrhip  of  the  planet -pass  in¬ 
duced  vibration  as  observed  at  a  fixed  point  on  'he  ring  gear.  The 
amplitude- frequency  spectrum  of  the  vibration  is  presented  by  Equations  (18) 
and  (19) . 


PLANET-PASS  INDUCED  VIBRATION  AND  SIDEBAND  AMPLITUDE  CALCULATIONS 


The  gear  parameters  for  the  CH-47  lower  planetary  gear  system  and  the  oper¬ 
ating  conditions  used  in  the  calculations  are  as  follows: 

Number  of  sun  gear  teeth  Ng  ■  28 

Number  of  planet  gear  teeth  N  ■  39 

P 


■> 


\ 
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Number  of  ring  gear  teeth  N^  ■  106 

Pitch  radius  of  sun  gear  r  *  2.8  in. 

s 

Pitch  radius  of  planet  gear  r  =  3.9  in. 

P 

Pitch  radius  of  ring  gear  rf  =  10.6  in. 

Pressure  angle  cp  **  25° 

Ring  gear  tooth  face  width  L  =  1.25  in. 

Angular  offset  of  planets  C  and  D  6  ■  1.343° 

Planet-ring  mesh  frequency  f^  *  1482  Hz 

Tangential  tooth  force  at  planet-ring  mesh  =  159.2  lb 

The  planet  orbiting  speed  is  14  Hz,  and  the  planet-pass  frequency  is 

therefore  56  Hz.  Due  to  the  angular  offset  6,  the  temporal  phase  lag  of 
planets  C  and  D  relative  to  A  and  B  is  37.6  degrees. 

The  cross  section  of  the  ring  gear  casing,  which  is  considered  as  the 

vibrating  elastic  body,  is  depicted  in  Figure  3.  In  the  dynamic  response 

calculation,  this  ring  gear  casing  is  modeled  as  a  composite  cylindrical 

shell  with  triable  thickness  (see  Figure  3).  Both  ends  of  the  casing  are 

assumed  to  be  "simply  supported,"  i.e.,  no  linear  translation  but  free  to 

rotate.  The  MTI  general  shell  dynamic  response  computer  program  was  used 

to  obtain  the  response  functions  b.  (9,z)  and  b_  (0,z)  [4]. 

AB  CD 

Let  the  vibration  observation  point  on  the  ring  gear  be  located  at  zq  =  1.5 

inches.  The  response  functions  at  this  point  for  the  specified  tooth  load 
are  plotted  in  Figure  4.  The  small  angular  offset  of  1.343°  for  planets  C 
and  D  is  neglected.  It  is  noted  that  both  b^g  and  bCD  are  periodic  with  a 

periodicity  of  180°.  Due  to  the  axial  symmetry  of  the  ring-gear  casing, 
b,  and  b  _  are  completely  similar  but  have  a  phase  difference  of  90°. 

The  vibration  induced  by  planet-pass  as  observed  at  a  fixed  point  located 
at  zq  -  1.5  inches  on  the  ring  gear  was  calculated  by  using  Equation  (17) 

and  is  plotted  in  Figure  5.  The  first  half  bump  corresponds  to  the  pass  of 
planet  A,  and  the  next  three  bumps  correspond  to  the  passes  of  planets  D,  B 
and  C,  respectively.  The  pattern  repeats  itself  after  one  full  revolution 
of  all  planets.  The  period  is  0.0714  sec,  which  is  the  reciprocal  of  u^. 

The  frequency  spectrum  of  this  planet-pass  induced  vibration  is  shown  in 
the  lower  diagram  of  Figure  6.  The  amplitude  at  the  mesh  frequency  has  the 
largest  value.  However,  the  amplitude  at  the  second  pair  of  sidebands,  at 
+  f  ,  is  quite  large. 

To  see  the  effect  of  planet  phasing  on  the  vibration  amplitudes  at  the 
selected  point,  the  vibrations  produced  at  three  different  temporal  phase 
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ACTUAL  CASING  CONFIGURATION 


Figure  3. 


Dynamic  Model  of  CH-47  Lower 
Ring-Gear  Casing. 


Planetary 


angles  (Y  *  0°,  90°  and  142.4°)  were  investigated.  The  amplitude-time  plots 
of  these  cases  are  shown  in  Figures  7,  8  and  9.  In  a  gross  sense,  they 
appear  to  be  similar.  However,  the  frequency  spectra  for  these  three  cases 
are  quite  different.  They  are  shown  in  the  upper  diagram  of  Figure  6  and 
the  upper  and  lower  diagrams  of  Figure  10  for  Y  *  0°,  90°  and  142.4°,  re¬ 
spectively.  It  is  noted  that  for  Y  =  0°,  sidebands  occur  only  at  frequencies 
of  £.  +  n  r  . 


gure  9.  Radial  Vibration  Amplitude  vs.  Time  at  Fixed  Locatii 
Forward  Rotor-Drive  Gearbox  Ring  Gear  for  t  =  142 .4 


ANALYSIS  OF  GEAR-MESH-INDUCED 
HIGH-FREQUENCY  VIBRATION  SPECTRA  AND  CALCULATIONS 


ANALYSIS  OF  GEAR  MESH  EXCITATION 


The  mesh  of  gtars  with  perfect  involute  profiles  can  induce  vibration, 
simply  because  of  tie  nonuniform  deflection  resulting  from  varying  tooth 
compliance  along  the  length  of  the  tooth.  The  nonuniform  deflection  or  dis¬ 
placement  deviation  introduces  irregular  motion  superimposed  on  the  uniform 
rotation  of  the  gears.  This  irregularity  of  motion  becomes  one  of  the  major 
sources  of  vibration  in  the  drive  system,  especially  of  the  vibration  asso¬ 
ciated  with  noise  production.  Manufacturing  errors  in  tooth  shape  or 
spacing  can  further  Increase  the  magnitude  of  the  motion  irregularity  and 
thus  the  excited  vibration  amplitude.  If  all  the  teeth  in  each  of  the 
meshing  gears  are  identical,  and  if  the  mesh  is  otherwise  ideal  (i.e.,  no 
runout,  etc.),  then  the  spectrum  of  the  vibration  induced  by  the  gear  mesh 
contains  signals  at  only  the  mesh  frequency  and  its  higher  harmonics.  The 
analysis  of  this  gear  mesh  excitation  was  performed  by  Laskin,  Orcutt  and 
Shipley  [ l]  . 

However,  if  the  gear  mesh  deviates  from  the  ideal  (e.g.,  gear  runout,  dy¬ 
namic  torque,  etc.),  or  if  there  is  deviation  in  tooth  profiles  as  the  gear 
mesh  continues  from  one  set  of  mating  teeth  to  the  other,  the  tooth  deflec¬ 
tion  pattern  will,  in  general,  vary  from  tooth  pair  to  tooth  pair.  The 
associated  vibration  will  then  contain  signal  components  at  other  than  the 
mesh  frequency  and  its  harmonics.  These  components  are  the  so-called  gear 
mesh  excitation  sidebands.  It  is  the  purpose  of  this  work  to  study  the 
gear-mesh-induced  vibration  spectra  considering  three  gear  parameters  as 
variables  between  me&n  cycles.  These  three  parameters  are: 

1.  Center-Line  Distance 


This  variation,  in  most  cases,  is  due  to  shaft  or  gear  runout.  The 
center- line  distance  is  treated  to  be  constant  within  a  tooth  mesh, 
but  it  can  vary  between  mesh  cycles.  This  simplifying  approximation 
is  possible  because  the  variation  in  center-line  distance  due  to  shaft 
or  gear  runout  is  slow  with  respect  to  tooth  mesh  frequency. 

2 .  Tooth  Load 

This  variation  may  be  caused,  for  instance,  by  the  dynamic  effects  of 
other  gear  meshes  in  the  drive  system.  The  tooth  load  can  also  be 
variable  between  calculation  points  within  one  mesh  cycle.  This 
accommodates  high-frequency  dynamic  forces. 

3.  Tooth  Support  Compliance 

This  is  the  compliance  in  addition  to  the  compliance  due  to  tooth 
bending,  shearing,  rotation,  and  contact  deformation  normally  existing 
during  gear  mesh.  It  may  be  produced  by  elastic  nonuniformity  in  the 
gear  structure  supporting  the  gear  teeth.  This  compliance  can  also 
be  variable  at  all  calculation  points  within  one  mesh  cycle. 

Each  of  these  three  variations  is  in  general  periodic  with  a  definite 
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frequency.  For  example,  the  center  distance  variation  caused  by  shaft  run¬ 
out  has  a  frequency  equal  to  the  rotational  speed  of  the  shaft.  In  each  of 
the  cases,  there  exists  a  base  period  for  the  tooth  oef lection  pattern 
which  is  equal  to  the  reciprocal  of  the  largest  common  factor  between  the 
mesh  frequency  and  the  frequency  (or  frequencies)  of  variation  of  the  gear 
parameter  (or  parameters).  During  one  base  period,  integer  multiples  of 
gear  mesh  cycles  and  of  parameter  variation  cycles  occur.  The  spectrum  of 
the  associated  vibration  is  thus  discrete,  with  the  fundamental  frequency 
equal  to  the  reciprocal  of  the  base  period.  The  frequency  spectrum  is 
obtained  by  a  Fourier  analysis  of  the  tooth  deflection  profile  over  one 
base  period. 

However,  if  the  variation  of  the  gear  parameter  is  random,  the  spectrum  of 
the  gear-mesh-induced  vibration  will  be  continuous.  It  may  be  obtained  by 
calculating  the  power  spectral  density  function  of  the  tooth  deflection 
data  over  a  long  period  of  time.  The  method  of  power  spectra  calculation 
is  included  in  Appendix  I. 

Within  one  mesh  cycle,  the  method  of  calculating  the  tooth  deflection  is 
directly  similar  to  that  in  [l]  .  The  limitations  and  assumptions  in  the 
calculation  are  listed  below: 

1.  Only  spur  gears  are  treated,  and  these  are  treated  only  for  the  two 
cases  of  (a)  an  external  gear  driving  an  external  gear  and  (b)  an 
external  geav  driving  an  internal  gear.  Straight  bevel  gears  may 
be  treated  in  an  approximate  manner  by  replacing  them  with  equiva¬ 
lent  spur  gears  by  Tregold's  Approximation  [  13]  . 

2.  The  working  portions  of  the  tooth  profiles  are  essentially  involute. 
Design  and  manufacturing  profile  deviations  are  small  enough  so  a.) 
not  to  affect  load  location,  load  direction,  or  tooth  stiffness. 

3.  There  is  no  tip  interference,  either  due  to  excessive  addendum 
length  or  due  to  tooth  deflection  under  load. 

4.  In  any  single  Interval  between  the  pitch  points  of  two  successive 
pairs  of  teeth,  contact  and  load  carrying  are  limited  to  the  two 
successive  pairs  of  teeth.  In  the  same  interval,  there  must  be  at; 
all  times  at  least  one  pair  of  teeth  in  contact  and  carrying  load. 
This  prevents  consideration  of  cases  where  the  contact  ratio  is 
less  than  one  or  more  than  two;  it  may  in  some  unusual  designs  also 
eliminate  cases  where  the  contact  ratio  has  certain  intermediate 
values . 

5.  The  load  is  assumed  to  be  transmitted  uniformly  across  the  face  of 
the  gear  except  for  normal  end  effects  in  stress  distribution. 

This  excludes  any  consideration  of  face  crowning,  lead  modifica¬ 
tion,  lead  manufacturing  error,  gear  windup,  or  nonuniform  deflec¬ 
tion  of  gear  supports. 

6.  All  variations  in  tooth  deflection  as  the  load  point  moves  along 
the  tooth  profile  either  are  confined  to  elastic  effects  on  the 
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tooth  alone  or  can  be  supplied  as  point-by-polnt  compliances  as 
part  of  the  input  data.  This  means  that  variations  such  as  might 
result  fron  the  deflection  of  thin  rims  are  not  calculated  directly 
by  the  analysis. 

7.  The  contact  deformation  is  assumed  to  be  independent  of  the  tooth 
surface  lubricating  film. 

The  above  analysis  is  incorporated  in  the  computer  program  GGEAR.  It  is 
obtained  by  modifying  and  extending  the  program  GEARO  reported  in  [l].  The 
overall  structure  of  the  modified  program  is  shown  in  Figure  11.  The  modi¬ 
fications  consist  of  the  creation  of  the  main  program  GGEAR  and  the  sub¬ 
routine  SPECT,  and  some  changes  in  subroutines  GEARO,  FOUR  and  PLT  to 
accommodate  the  extension  of  computation  over  multiple  mesh  cycles.  GGEAR 
accepts  those  items  of  gear  data  that  are  constant  over  all  mesh  cycles. 

Subroutine  GEARO,  on  the  other  hand,  reads  in  variable  gear  data  and,  to¬ 
gether  with  subroutines  AJCDH  and  CALCJ,  calculates  the  tooth  deflection 
over  one  mesh  cycle.  In  the  computation,  the  mesh  cycle  is  divided  into  a 
number  of  calculation  points  as  is  done  in  the  program  GEARO.  Tooth  deflec¬ 
tions  at  all  the  calculation  points  over  the  prescribed  number  of  mesh 
cycles  are  stored  and  printed  out  in  GGEAR. 

With  reference  to  Figure  11,  it  may  be  noted  that  provision  is  being  made 
for  eventual  incorporation  o£  other  subroutines  similar  to  GEARO  for  the 
calculation  of  excitation  levels  in  high-contact-ratio  Bpur  gears,  helical 
gears,  and  eventually  spiral-bevel  gears.  While  such  calculation  capabili¬ 
ties  were  not  included  in  the  computer  program  during  the  present  contract, 
the  modified  program  was  prepared  on  a  modular  basis  for  their  later  inclu¬ 
sion. 

According  to  the  user's  instruction,  either  the  calculated  tangential  deflec¬ 
tion  data  can  be  plotted  by  subroutine  PLT,  or  it  can  be  analyzed  to  obtain 
the  Fourier  representation  of  the  deflection  pattern  via  the  subroutine  FOUR, 
as  is  now  done  in  program  GEARO.  The  data  can  also  be  analyzed  by  using 
SPECT  to  produce  the  power  spectral  density  function.  A  description  of  the 
computer  program  is  given  in  Appendix  II . 

CH-47  SPIRAL  BEVEL  GEAR  MESH  CALCULATIONS 

The  CH-47  forward  rotor-drive  gearbox  spiral  bevel  gear  mesh  is  analyzed  in  *• 

terms  of  its  equivalent  spur  gear  mesh  [2].  The  equivalent  spur  gears  will  JW 

be  equivalent  to  the  spiral  bevel  gears  only  in  the  sense  that  their  . 

physical  proportions,  those  likely  to  influence  deflection  under  load,  ; 

approximate  the  mean  proportions  of  the  actual  spiral  bevel  gear  teeth.  The 
use  of  the  equivalent  spur  gears  is  sufficiently  adequate  for  the  present 
purpose  of  studying  the  effects  of  shaft  runout  and  load  variation  to  pro¬ 
duce  high- frequency-vibration  sidebands.  The  gear  mesh  parameters  and  the 
conversion  of  the  CH-47  spiral  bevel  gears  into  equivalent  spur  gears  are 
summarized  in  Table  I,  which  has  been  taken  from  [2].  The  gear  mesh  'i 

tangential  load  is  taken  to  be  2760  lb,  and  the  gearbox  input  shaft  speed  is 
7059  rpm.  This  yields  a  gear  tooth  mesh  frequency  of  3412  Hz.  Variations 
in  center  distance  and  in  tooth  load  as  the  source  for  generating  tooth  dis¬ 
placement  dev.  .itions  at  sideband  frequencies  have  been  studied. 
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Figure  11.  Computer  Program  Structure 


1.  Center  Distance  Variation 


Center  distance  variation  resulting  from  a  0.001-in.  (peak)  input  shaft 
runout  is  chosen  to  be  the  representative  case.  The  runout  is  assumed 
to  be  a  siie  wave  in  form.  The  expression  for  the  center  distance  is 

2nf 

C  (in.)  -  22.3735  +  0.001  sin  t  (21) 

where  the  number  22.3735  is  the  nominal  center  distance  in  inches,  t  the 
time  in  seconds,  and  f^  the  mesh  frequency  in  Hz.  The  input  pinion  has 

29  teeth,  and  therefore  f^/29  is  the  shaft  rotational  speed. 

Using  the  computer  program,  tooth  displacement  deviation  is  calculated 
over  29  mesh  cycles  covering  one  period  of  variation  in  center  distance 
(one  complete  shaft  rotation).  The  tooth  deviation  is  periodic  with  a 
period  equal  to  one  shaft  rotation.  This  periodic  tooth  deviation  pat¬ 
tern  is  analyzed  by  using  the  extended  subroutine  FOUR  to  obtain  the 
amplitude- frequency  relationship.  This  relationship  is  shown  in  Figure 
12.  Amplitude  is  the  magnitude  of  tangential  deviation  over  and  above 
the  mean  tooth  deflection.  The  frequency  is  expressed  in  multiples  of 
mesh  frequency  f^.  Since  the  tooth  deviation  is  periodic,  the  amplitude 

distribution  is  discrete  and  is  nonzero  only  at  the  multiples  of  re¬ 
ciprocal  of  the  base  period,  which  is  the  shaft  speed  (f^/29).  It  is 

seen  from  Figure  12  that  mesh  frequency  is  still  the  most  dominant 
frequency.  However,  the  sideband  amplitudes  close  to  the  harmonic  fre¬ 
quencies  are  seen  to  be  comparable  to  the  high  harmonic  components. 

It  is  expected  that  a  larger  shaft  runout  would  produce  higher  sideband 
amplitudes.  This  is  indeed  shown  in  Figure  13,  which  is  a  frequency- 
amplitude  plot  for  a  0.002-in.  (peak)  shaft  runout.  It  is  seen  that  the 
amplitudes  at  f^/29  and  at  sidebands  around  the  higher  harmonics  are 

approximately  doubled  in  magnitude. 

2.  Transmitted  Load  Variation 


A  sinusoidal  dynamic  load  of  500  lb  at  half  mesh  frequency  is  assumed 

to  be  superimposed  on  the  nominal  tooth  load  of  2760  lb.  The  total  » 

tangential  tooth  load  is  therefore 

2rf>‘  * 
W  (lb)  -  2760  +  500  sin  t  (22) 

The  calculated  tooth  displacement  deviation  has  a  periodicity  of  two 

mesh  cycles.  The  resulting  amplitude- frequency  plot  of  this  periodic 

tooth  deviation  pattern  is  shown  in  Figure  14.  It  is  seen  that  the 

first  sideband  amplitude,  which  is  at  fM/2,  is  larger  than  all  the  high 

**  » 

harmonic  components.  Also,  all  other  sideband  amplitudes  are  relatively 

small.  This  large  first  sideband  amplitude  may  be  attributed  to  the 

large  maximum  dynamic  load  (500  lb)  relative  to  the  nominal  load. 
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CH-47  Spiral  Bevel  Gear  Mesh. 


Figure  14.  Frequency  vs.  Excitation  Amplitude  for  Tooth  Force 
Variation,  CH-47  Spiral  Bevel  Gear  Mesh. 


CH-47  LOWER  PLANETARY  PLANET-TO-RING  GEAR  MESH  CALCULATIONS 


The  CH-47  lower  planetary  planet  and  ring  gear  parameters  are  summarized  in 
Table  II.  The  tangential  tooth  load  at  the  planet-to-ring  gear  mesh  is 
taken  to  be  2103  lb.  The  planet  rotational  speed  is  2280  rpm.  The  corre¬ 
sponding  mesh  frequency  is  1482  Hz.  Effects  of  variation  in  center 
distances,  tooth  load,  and  tooth  support  compliance  to  produce  vibration 
sideband  frequencies  have  been  investigated. 

1.  Center  Distance  Variation 


A  planet  runout  of  0.001  in.  (peak)  in  sine  wave  form  is  assumed.  With 
a  nominal  center  distance  of  6.7  in.,  the  expression  for  center  distance 
is 

2TTfM 

C  (in.)  =  6.7 +0.001  sin  —■  t  (23) 

Tooth  displacement  deviation  is  calculated  using  the  computer  program, 
over  39  mesh  cycles,  to  cover  a  full  rotation  of  the  planet.  The  fre¬ 
quency  spectrum  of  the  tooth-mesh -induced  vibration,  obtained  by 
Fourier  analyzing  the  displacement  deviation,  is  plotted  in  Figure  15. 
The  mesh  frequency  is  still  the  dominating  frequency.  The  largest 
sideband  amplitudes  occur  at  the  sidebands  closest  to  and  on  each  side 
of  the  mesh  frequency. 


TABLE  II.  CH-47  LOWER  PLANETARY  PLANET  AND 

RING  GEAR  PARAMETERS 

Planet 

Ring 

Number  of  Teeth 

39 

106 

Face  Width 

1.55  in. 

1.25 

in. 

Pitch  Radius 

3.9  in. 

10.6 

in. 

Outer  Radius 

4.0845  in. 

- 

Inner  Radius 

- 

10.43 

in. 

Root  Radius 

3.69  in. 

10.845 

in. 

Radius  to  the  Beginning  of  Involute  Profile 

3.738  in. 

10.7935 

in . 

Circular  Tooth  Thickness  at  Pitch  Circle 

0.3462  in. 

0.276 

in. 

Tooth  Fillet  Radius 

0.075  in. 

0.094 

in. 

Pressure  Angle 

25  deg 

25  deg 

i  . 

Figure  15.  Frequency  vs.  Excitation  Amplitude  for  Planet  Runout,  CH-47 
Lower  Planetary  Planet-to-Ring  Gear  Mesh. 


2  .  Transmitted  Load  Variation 


A  dynamic  load  of  400  lb  varying  sinusoidally  at  one-third  of  the  mesh 
frequency  is  assumed  to  be  superimposed  on  the  nominal  transmitted  load  of 
2103  lb.  The  total  tooth  load  is  therefore 

2lTfM 

W  (lb)  =  2103  +  400  sin  t  (24) 

The  vibration  induced  by  this  load  variation  at  gear  tooth  mesh  has  a 
fundamental  frequency  of  f^/3.  The  frequency  spectrum  is  shown  in 

Figure  16,  which  is  obtained  by  the  computer  program  calculating  the 
tooth  displacement  deviation  over  three  mesh  cycles.  It  is  seen  that 
the  amplitude  at  f^/2  is  higher  than  those  at  high  harmonics  of  the 
mesh  frequency. 

3.  Tooth  Support  Compliance  Variations 

Five  consecutive  ring  gear  teeth  (out  of  a  total  of  106  teeth)  are 
assumed  to  have  a  finite  (tangential)  tooth  support  compliance  of  10"° 
in. /lb.  This  compliance  is  in  addition  to  the  tooth  bending,  tooth 
shear,  tooth  rotation,  and  contact  deformation  compliances  calculated  by 
the  computer  program  in  terms  of  gear  tooth  geometry  and  elastic  prop¬ 
erties  . 

The  amplitude- frequency  plot  of  the  vibration  induced  by  this  compliance 
variation  is  shown  in  Figure  17.  The  information  shown  in  this  figure 
may  be  Interpreted  as  follows:  Assume  first  that  the  planets  are  fixed 
and  the  ring  gear  rotates.  The  frequency  of  mesh  between  a  fixed  planet 
and  the  ring  gear  is  equal  to  the  number  of  teeth  on  the  ring  gear 
multiplied  by  the  relative  rotational  speed  between  them.  This  speed  is 
actually  the  orbiting  speed  of  the  planet  with  respect  to  a  fixed  body, 
since  the  ring  gear  is  stationary.  The  periodicity  of  the  tooth  deflec¬ 
tion  excitation  caused  by  an  assumed  local  ring  gear  compliance  varia¬ 
tion  is  equal  to  the  time  for  one  full  relative  rotation  of  the  ring 
gear  with  respect  to  the  planet.  Since  there  are  106  teeth  on  the  ring 
gear,  this  periodicity  is  106  times  the  tooth  mesh  time;  that  is, the 
reciprocal  of  the  tooth  mesh  frequency.  Therefore,  in  the  amplitude- 
frequency  plot  shown  in  Figure  17,  there  appear  105  equally-spaced  side¬ 
band  frequencies  between  two  successive  harmonics  of  mest  frequency. 

Any  two  neighboring  sidebands  are  separated  by  the  planet  orbiting  fre¬ 
quency  . 

It  is  seen  from  the  figure  that  the  sideband  amplitudes  are  in  general 
quite  large  and  that  the  amplitudes  of  the  low-frequency  sidebands  are 
even  larger  than  the  amplitude  at  the  mesh  frequency.  The  compliance 
due  to  tooth  bending,  tooth  shear,  tooth  rotation  and  contact  deforma¬ 
tion  is  about  10  in. /lb  calculated  by  the  computer  program.  The  large 
sideband  amplitudes  are  therefore  due  to  the  large  tooth  support  com¬ 
pliance  variation  relative  to  the  normal  compliance  existing  during 
tooth  mesh. 
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EXCITATION  AMPLITUDE  TANGENT  TO  PITCH  CIRCLE  (I0‘*  IN.) 
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It  is  also  seen  from  Figure  17  that  sidebands  are  grouped  into  a  series 
of  bumps,  and  there  are  about  five  bumps  between  two  adjacent  harmonic 
frequencies.  This  is  attributed  to  the  fact  that  the  finite  tooth  sup¬ 
port  compliance  variation  occurs  in  only  five  consecutive  teeth  and 
that  the  compliance  variation  is  much  larger  than  the  normal  tooth  mesh 
compliance.  Tooth  deflection  patterns  over  one  period  (106  mesh  cycles) 
may  be  considered  as  the  sum  of  two  separate  deflection  profiles.  One 
is  the  normal  tooth  mesh  deflection  profile  consisting  of  106  identical 
segments  representing  deflections  for  the  106  mesh  cycles.  The  other 
is  the  deflection  due  to  the  tooth  support  compliance  variation,  which 
is  nonzero  only  in  the  first  five  mesh  cycles  and  is  much  larger  in 
magnitude  than  the  former.  The  frequency  content  of  the  first  deflec¬ 
tion  profile,  consists  of  only  the  mesh  frequency  and  its  harmonics, 
while  the  frequency  spectrum  of  the  latter  is  roughly  periodic  with 
declining  ampliti.de  and  the  period  is  about  one-fifth  of  the  mesh  fre¬ 
quency.  Therefore,  the  frequency  spectrum  of  the  total  tooth  deflection 
pattern  appears  to  be  a  series  of  bumps  superimposed  on  the  mesh  fre¬ 
quency  harmonic 8. 

CH-47  SPIRAL  BEVEL  GEAR  SYSTEM  RESPONSE  CALCULATIONS 


In  order  to  demonstrate  the  manner  in  which  the  gear  excitation  sidebands 
may  be  applied  in  the  vibration  analysis  of  geared  power  trains,  and  to 
illustrate  the  type  of  response  which  is  produced  by  gear  excitation  spec¬ 
tra  which  include  sidebands,  a  sample  calculation  has  been  made  using  com¬ 
puter  program  TORRP  [2] . 

In  this  calculation,  the  excitation  spectrum  obtained  under  the  assumption 
of  0.002-in.  runout  in  the  CH-47  spiral  bevel  mesh  was  used.  Mesh  fre¬ 
quency  is  3412  Hz,  and  a  steady-state  tangential  tooth  force  of  2760  lb 
was  used.  The  gearbox  drive  train  components  were  represented  dynamically 
by  the  torsional  response  model  reported  in  [2]  . 

Results  of  the  calculations  are  shown  in  Table  III  and  Figures  18  and  19. 
From  these  results  it  is  apparent  that  significant  torsional  response  occurs 
in  the  gearbox  components  as  a  result  of  the  sideband  disturbances.  It  is 
worth  mentioning  that  essentially  zero  response  occurs  at  frequencies 
between  the  peaks  shown  because  there  are  no  sources  of  excitation  at  these 
frequencies . 


TABLE  III  -  Continued 


.002-IN.  SHAFT  RUfKX 
CH-47  SPIRAL.  BtVEL 


Figure  18.  Calculated  Dynamic  Tooth  Force  at  CH-47  Spiral  Bevel  Gear  Mesh 


Figure  19.  Calculated  Dynamic  Tooth  Force  at  CH-47  Planet-Ring  Gear  Mesh. 


DISCUSSION  OF  RESULTS 


DISCUSSION  OF  TASK  I  RESULTS 


The  planet-pass  Induced  vibration  given  by  Equation  (17)  is  the  amplitude¬ 
time  history  of  radial  displacement  at  a  fixed  point  on  the  ring  gear 
casing  The  zero  time  is  arbitrarily  chosen  to  be  the  moment  when  planet 
A  (see  Figure  2)  is  at  the  observation  point.  As  time  progresses,  planets 
D,  B  and  C  pass  the  observation  point  in  sequence.  This  sequence  repeats 
with  a  frequency  equal  to  the  planet  orbiting  speed  a^. 

Since  ring  gear  radial  vibration  results  from  the  planet-ring  gear  mesh 
dynamic  tooth  force,  the  vibration  amplitude  is  large  when  the  planet  is 
close  to  the  observation  point  and  smaller  when  the  planet  is  more  distant. 
This  is  seen  from  the  calculated  results  in  Figures  5,  7,8  and  9,  which 
are  for  phase  angles  i|t  =*  37.6°,  OP,  90°  and  142.4°  (between  opposite  pairs 
of  planets),  respectively.  The  four  major  large-vibration-amplitude  areas 
in  each  of  the  figures  correspond  to  the  passage  of  the  four  planets  past 
the  fixed  observation  point. 


In  general,  it  appears  that  the  shape  and  magnitude  of  the  vibration  data 

in  Figures  5  through  9  are  quite  similar,  and  that  only  the  amplitude 

between  planet  passes  varies  appreciably  with  the  phase  angle  t|i .  This  is 

because  the  two  components  of  the  total  response  function,  b  and  b  , 

^  A8  CD 

are  quite  localized  (i.e.,  peaked  at  0  and  90°,  respectively).  Also, 

they  are  identical  in  shape  due  to  the  axisymmetry  of  the  ring  gear  casing, 

but  with  90°  phase  difference  (see  Figure  4).  This  90°  phase  difference 

corresponds  to  the  time  spacing  between  the  passes  of  adjacent  planets. 

Us.ng  the  first  two  bumps  as  examples,  the  vibration  amplitude  around  the 

passes  of  planets  A  and  D  (around  9  =  0°  end  90°)  are  dominantly  determined 

by  b^B  and  b^,  respectively,  without  much  interference  between  them. 

Only  between  planet  passes  (around  9  =  45°)  does  appreciable  interference  exist, 
and  therefore  the  local  amplitude  depends  on  the  temporal  phase  difference 
i|r  between  the  two  response  functions  (see  Equations  (6)  and  (7)). 


The  vibration  of  the  ring  gear  casing  is  induced  by  the  passing  of  the 
planet  gears  and  originates  in  the  dynamic  tooth  forces  which  exist  in 
the  planet-ring  gear  tooth  meshes.  The  vibration  occurs,  therefore,  at 
the  gear  tooth  mesh  frequency,  f^,  modulated  by  the  planet-pass  frequency. 


This  is  shown  in  the  frequency  spectrum  plots  of  Figures  6  and  10.  It  is 
seen  that  the  mesh  frequency  is  the  center  frequency  and  that  around  it, 
there  are  a  number  of  sidebands  which  are  found  at  mesh  frequency  plus  or 
minus  (in  general)  integer  multiples  of  half  of  the  planet-pass  frequency 
fp.  In  the  special  case  of  if  ■  0,  four  equally  spaced  planets  (small 

positional  offset  6  has  been  neglected)  pass  the  observation  point  with 
orbiting  speed  ui  .  Dynamically,  this  is  equivalent  to  one  planet  passing 


this  point  wieh  planet-pass  frequency  f^  (where  f^  ■  4  o^) .  Therefore, 
sidebands  occur  only  at  f^  ±  n  f  ,  where  n  is  an  integer.  In  the  more 


general  case  where  \|i  t  0,  there  is  some  phase  difference  between  the 
planet  pairs  (A,B)  and  (C,D).  Each  pair  passes  the  observation  point 
with  a  frequency  of  f^/2.  Therefore,  sidebands  in  general  appear  at 

£m  *  n  f  ■ 

Comparison  of  the  four  spectrum  diagrams  shown  in  Figures  6  and  10  reveals 
that  the  effect  of  the  phase  angle t  is  to  increase  the  sideband  amplitude 
at  mesh  frequency  plus  and  minus  integer  multiples  of  f^/2.  For 

i|r  =  0,  the  amplitudes  at  thesa  sidebands  are  zero,  while  for  i(r  =  142.4°, 
the  corresponding  amplitudes  are  even  greater  than  those  at  the  neigh¬ 
boring  frequencies  of  ^  i  n  f  . 

Under  some  circumstances,  it  may  be  desirable  to  control  the  vibration 
amplitudes  at  certain  sidebands.  For  the  CH-47  lower  planetary  gear 
system,  the  phase  angle  ^  is  a  major  controlling  parameter  of  sideband 
amplitudes  for  the  planet-pass  induced  vibration.  The  angle  in  turn  is 
determined  by  other  gear  parameters  including  the  numbers  of  teeth  of 
component  gears.  For  example.  Equation  (4)  indicates  that  the  greater 
the  number  of  teeth  on  the  sun  gear,  the  larger  the  phase  angle.  However, 
other  gear  performance  characteristics  such  as  speed  ratio,  load  distri¬ 
bution,  etc.,  must  be  considered  when  gear  parameters  are  altered  to  con¬ 
trol  the  value  of  ^ - 

Finally,  it  should  be  noted  that  the  analysis  used  herein  assumes  that  the 
normal  components  of  the  planet  gear  dynamic  tooth  forces  vary  sinusoid¬ 
ally  with  time.  This  variation  occurs  about  the  nominal  steady-state 
normal  tooth  force  component  which  accompanies  the  transfer  of  torque. 
Passage  of  this  nominal  steady-state  force  itself  causes  quasi-static 
ring  gear  deflections,  which  are  ignored  in  the  vibration  calculations. 

DISCUSSION  OF  TASK  II  RESULTS 


Vibration  can  be  excited  by  nonuniform  tooth  deflection  and  tooth  pro¬ 
file  characteristics  during  the  mesh  of  even  precisely  machined  gears. 

If  the  center  distance  and  nominal  tooth  load  are  constant  and  all  other 
gear  parameters  are  the  same  from  one  pair  of  teeth  to  the  next,  the 
nonuniform  tooth  deflection  pattern  in  one  mesh  cycle  will  repeat  in  all 
successive  cycles.  The  frequency  content  of  this  tooth-mesh- induced  vi¬ 
bration  will  consist,  therefore,  of  only  the  mesh  frequency  and  its  har¬ 
monics.  However,  if  there  is  any  variation  in  center  distance  or  tooth 
load,  or  any  change  of  any  other  gear  parameters  from  mesh  to  mesh,  other 
vibration  frequencies  (the  so-called  sidebands)  are  introduced.  The 
results  of  this  investigation  have  clearly  shown  this  point. 

From  the  frequency  spectra  shown  in  Figures  12,  13,  and  15,  it  is  seen 
that  major  sidebands  due  to  variation  in  center  distance  occur  around  the 
mesh  frequency  and  its  harmonics.  In  all  cases  the  mesh  frequency  is 
still  the  most  dominating  frequency.  A  comparison  of  Figures  12  and 
13  reveals  that  the  sideband  amplitudes  depend  quite  strongly  on  the 
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magnitude  of  the  center  distance  variation.  The  general  shape  and  magni¬ 
tude  of  the  harmonic  components  in  Figure  13  (0.002-in.  runout)  are  similar 
to  those  shown  in  Figure  12  (0.001-in.  runout).  However,  the  amplitudes 
at  major  sidebands  in  Figure  13  are  about  twice  those  at  corresponding 
sidebands  shown  in  Figure  12.  Furthermore,  in  the  0.002-in.  runout  case, 
the  sideband  amplitudes  around  the  fourth  and  fifth  harmonics  are  even 
larger  than  the  harmonic  amplitudes. 

Actual  vibration  frequency  spectrum  measurements  [8]  taken  from  an 
accelerometer  mounted  on  the  CH-47  rotor-drive  gearbox  bear  some  similar¬ 
ity  to  those  in  Figures  12,  13,  and  15.  Output  from  this  accelerometer, 
located  on  the  outside  of  the  gearbox  near  the  spiral  bevel  gear  shaft 
support  bearings,  is  reproduced  in  Figure  20.  It  should  be  noted  that  the 
noise  signals  shown  in  Figure  20  are  generated  not  only  by  the  spira? 
bevel  gear  mesh  but  also  by  a  number  of  other  possible  vibration  sources, 
such  as  gear  runout,  planet-pass,  etc.,  occurring  in  the  drive  system. 

The  signal  is  strongest  at  the  spiral  bevel  mesh  frequency  because  the 
accelerometer  is  located  near  the  spiral  bevel  gears' shaft  support  bear¬ 
ings.  The  sideband  with  large  amplitude  near  the  spiral  bevel  gear  mesh 
fundamental  is  separated  by  118  Hz  from  the  mesh  frequency.  This  differ¬ 
ence  is  the  rotational  speed  of  the  input  pinion.  Therefore,  this  side¬ 
band  is  clearly  due  to  input  shaft  runout, and  it  is  equivalent  to  the 
relatively  large-amplitude  sidebands  around  the  mesh  frequency  shown  in 
Figures  12  and  13. 

The  sideband  amplitudes  produced  by  tooth  load  variations  as  shown  in 
Figures  J4  and  16  are  quite  large  in  comparison  with  the  harmonic  compo¬ 
nents.  In  both  cases  the  amplitude  of  the  first  sideband  is  greater  than 
the  amplitudes  at  high  harmonic*  of  the  mesh  frequency.  These  large 
sideband  amplitudes  may  be  attributed  to  the  large  load  variation  relative 
to  the  nominal  tooth  load  (500  lb  vs.  2760  lb  for  spiral  bevel  gear  mesh 
and  400  lb  vs.  2103  lb  for  lower  planetary  planet-to-ring  gear  mesh). 

In  the  example  of  tooth  support  compliance  variatijn,  five  consecutive 
teeth  were  assumed  to  have  a  finite  compliance  of  10"  in. /lb  in  the 
tangential  direction.  This  represents  a  case  in  which  the  lower  planetary 
ring  gear  casing  has  a  relatively  weak  local  elastic  stiffness.  This 
finite  support  compliance  produces  tangential  tooth  displacement  in  addi¬ 
tion  to  that  due  to  normal  tooth  bending,  shear,  rotation  and  contact  de¬ 
formation  during  mesh.  Figure  17  shows  that  the  sideband  amplitudes  re¬ 
sulting  from  this  compliance  consist  of  a  number  of  bumps  superimposed  on 
the  harmonic  components.  The  amplitudes  at  small  sideband  frequencies  are 
large, and  those  at  the  first  few  sidebands  are  even  greater  than  the  ampli¬ 
tude  at  the  gear  mesh  frequency.  This  is  because  the  nominal  compliance 
with  respect  to  tooth  deflection  during  gear  mesh  is  only  on  the  order  of 

10  ^  in. /lb,  which  is  one  order  of  magnitude  smaller  than  the  assumed 
local  tooth  support  compliance. 

In  all  of  the  cases  studied  in  this  investigation,  the  frequency  spectrum 
is  discrete  (i.e.,  the  amplitudes  are  nonzero  only  at  discrete 
frequencies),  because  the  frequency  of  variation  of  the  gear 
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parameter  selected  for  study  Is  In  all  cases  a  rational  fraction  of  the 
gear  mesh  frequency.  There  thus  always  exists  a  base  period  over  which 
the  tooth  displacement  pattern  repeats  itself.  This  base  period  is  an 
integer  multiple  of  the  gear  mesh  time,  which  is  the  reciprocal  of  the 
mesh  frequency.  Therefore,  the  vibration  frequencies  are  integer  multi¬ 
ples  of  the  reciprocal  of  the  base  period, and  the  frequency  spectrum  is 
thus  discrete.  The  frequency  spectrum  is  obtained  by  performing  a  Fourier 
analysis  of  the  periodic  tooth  deflection  pattern. 

If  the  variation  of  the  gear  parameter  is  random,  or  if  its  frequency  is 
an  irrational  fraction  of  the  mesh  frequency,  the  corresponding  frequency 
spectrum  of  the  gear-mesh- induced  vibration  will  be  continuous.  The  sub¬ 
routine  SPECT  may  be  used  to  compute  the  power  spectral  density  function 
based  on  the  tooth  deflection  data  over  a  large  number  of  mesh  cycles  in 
such  cases. 
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CONCLUSIONS 


The  following  conclusions  are  drawn  as  a  result  of  the  studies  reported 
herein: 

With  respect  to  planet-pass  vibration  sidebands: 

1.  Methods  have  been  developed  for  predicting  geometrically- induced 
planet-pass  vibration  sidebands  which  accompany  normal  planetary 
gear  reduction  operation. 

2.  Planet-pass  vibration  sideband  amplitudes  may  exceed  that  of  the 
base  signal,  and  plans  for  dealing  with  them  must  be  part  of  any 
vibration  and  noise  reduction  program. 

3.  Planet-pass  vibration  sideband  frequencies  exist  both  below  and 
above  the  base  signal  at  integer  multiples  of  one-half  planet- 
pass  frequency  in  four-planet  systems  with  opposite  pairs  of 
planets  in  phase. 

4.  Planet-pass  vibration  sideband  spectra  are  affected  by  the  phase 
relationships  of  the  planets,  and  spectra  for  several  such  rela¬ 
tionships  have  been  determined  for  the  CH-47  forward  rotor  drive 
gearbox  first-stage  planetary  reduction. 

With  respect  to  vibration  sidebands  in  single  gear  meshes: 

1.  Methods  have  been  developed  for  predicting  vibration  sidebands 
which  are  produced  by  gear  runout,  dynamic  variations  in  tooth 
transmitted  force,  and  tooth  support  discontinuities. 

2.  Shaft  runout  vibration  sideband  amplitudes  may  exist  with  signif¬ 
icant  amplitudes.  They  must  therefore  be  considered  in  vibration 
and  noise  reduction  efforts.  Conversely,  such  sidebands  if  detect¬ 
able  could  be  used  for  the  diagnosis  of  improperly  assembled  shafts 
and  bearings. 

3.  Dynamic  transmitted  tooth  force  variations  can  produce  vibration 
sidebands  with  significant  amplitudes.  While  the  presence  of  such 

signals  is  important  from  the  vibration  and  noise  reduction  stand-  ( 

point,  the  most  important  aspects  of  this  type  of  sideband  are  as  ^ 

follows:  * 

a)  This  type  of  disturbance  implies  the  presence  of  additional 
dynamic  loads  on  gear  teeth  and  bearings,  with  resulting 
degradation  of  component  lifetimes.  The  importance  of  this 

type  of  dynamic  force  on  gear  and  bearing  life  has  yet  to  be  ,'v 

assessed . 

b)  A  major  source  of  vibration  signals  within  operating  gearboxes 
has  been  explained,  thereby  reducing  the  number  of  unknown 
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signals  which  must  be  explained  in  any  diagnostic  exercise. 
Likewise,  the  dynamic  interactions  between  gear  meshes  in 
multiple-mesh  gear  trains  have  been  explained. 

4.  Vibration  sidebands  of  a  very  distinctive  character  are  predicted 
herein  to  be  produced  by  gear  meshes  with  tooth  support  discontin¬ 
uities,  such  as  cracks.  Consequently,  a  very  important  diagnostic 
method  has  been  identified  for  identifying  gear  structural  in¬ 
tegrity  problems  which  may  be  introduced  during  manufacture  or 
which  may  appear  in  service. 

In  summary,  methods  have  been  devised  both  for  designing  low-vibration  and 
noise  gear  reductions  and  for  identifying  the  existence  of  several  types 
of  gear  problems.  It  is  felt  that  the  diagnostic  potential  of  the  vibra¬ 
tion  analysis  methods  described  herein  is  of  considerable  importance. 


RECOMMENDATIONS 


The  results  reported  herein  strengthen  considerably  an  already  valuable 
body  of  technology  which  relates  gear  train  component  condition  to  measur¬ 
able  symptoms.  The  most  important  aspect  of  this  technology  is  that  much 
of  it  appears  as  computer-implemented  analytical  procedures  which  can  be 
utilized  by  the  gearbox  designer. 

Recommendations  for  future  efforts  in  this  important  area  fall  into  three 
specific  areas: 

1.  It  is  recommended  that  the  analytical  procedures  next  be  extended 
to  include  high-contact-ratio  gearing  and  coupled  response  of 
gear  train  components. 

2.  It  is  recommended  that  the  high-frequency  vibration  analysis  tools 
be  utilized  in  their  present  form,  where  applicable,  in  the  design 
of  future  gearboxes;  further,  that  the  results  of  this  usage  be 
documented  in  specific  instances,  particularly  where  comparative 
testing  accompanies  this  usage,  for  feedback  to  this  technology 
program. 

3.  It  is  recommended  that  the  use  of  high-frequency  gearbox  vibra¬ 
tion  technology  be  included  as  an  integral  part  of  future  efforts 
directed  at  drive  train  condition  monitoring  and  diagnosis,  and  at 
component  life  and/or  failure  prognosis. 
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APPENDIX  I 

CALCULATION  OF  POWER  SPECTRAL  DENSITY  FUNCTION 


The  method  of  calculation  presented  herein  follows  that  given  in  [12].  Let 
n  be  the  number  of  calculation  points  in  one  tooth  mesh  cycle.  The  time 
interval  between  data  points  is 


At 


(25) 


where  £,  is  the  mesh  frequency.  The  cutoff  frequency  is  defined  as 


f  * 

c 


(26) 


The  total  time  of  record  is  defined  as 


Tr  “  *  n  *  At 


(27) 


where  n^  is  the  number  of  mesh  cycles  over  which  tooth  deflection  data  is 
recorded. 


The  autocorrelation  function  is,  by  definition, 

T 

I  f 


R(t)  «*  j&im  —  ^  x  (t)  x  (t  +  T)  dt 


(28) 
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where  x  is  the  tooth  deflection  and  t  la  the  time.  The  power  spectral  density 
function  is  related  to  R(t)  by 


G  (f)  ■  4  \  R(t)  cos  2tt  fT  dT 


(29) 


where  f  is  frequency. 

Tooth  deflection  x  is  known  numerically  at  discrete  points  separated  by  At. 
Let 


y 


\ 


»  x(i  At)  i  =»  1,  2,  ...,  (30) 

The  autocorrelation  function  at  displacement  r  At  may  be  estimated  by 
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where  r  Is  the  lag  number  and  m  the  maximum  lag  number.  This  maximum  lag 
number  determines  the  maximum  displacement  and  the  equivalent  resolution 
bandwidth  for  power  spectral  calculation  as  follows: 


max 


m  At 


max 


m 

n  ^ 
Uh 
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(32) 


(33) 


It  is  desirable  to  keep  T  less  than  one-tenth  of  the  time  of  record  T  . 

max  r 

This  will  avoid  certain  instabilities  that  can  occur  in  autocorrelation 

function  estimates.  In  the  calculation. T  is  set  to  be  about  1/20  of  T  , 

*  max  r 

m  is  then  determined  by 

n  i\, 

(34) 
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In  other  words,  m  is  set  equal  to  the  integer  part  of  the  right-hand  side. 
The  numerical  approximation  of  Equation  (29)  is 

m-1 

^  f  f  I  _ .  I  nn  fkt 

(35) 
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The  numerical  estimate  of  the  power  spectral  function  should  be  calculated 
only  at  the  nrt-1  special  discrete  frequencies  where 


k  f 


k  *  0,  1,  2,  ...  m 


(36) 


This  will  provide  m/2  independent  spectral  estimates  since  the  bandwidth 
2  f 

B  is  ■-  ■  .  At  these  discrete  frequencies, 
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The  index  k  is  called  the  harmonic  number.  G,  is  the  "raw"  estimate  of  the 

k 

power  spectral  density  function  at  harmonic  k.  The  "smooth"  estimate  at 
harmonic  k  is 

G  -  0.5  G  +  0.5  G. 

o  o  1 

Gk  =  0.25  Gk-1  +  0.5  Gk  +  0.25  Gfcfl  k  =  1,  2,  ...m-1  (38) 

G  -  0.5  G  ,  +  0.5  G 

m  m-1  m 


The  ratio  of  the  frequency  at  harmonic  k  to  mesh  frequency  is 

£k  kn 
fM  =  2m 


(39) 


where  m  is  given  by  Equation  (34).  The  frequency  interval  is  approximately 


(40) 


The  above  method  of  calculation  has  been  progranmed  in  subroutine  SFECT. 
This  subroutine  has  been  tested  successfully  by  using  a  simple  sine  func¬ 
tion  over  a  hundred  cycles  of  time.  A  sharp  peak  was  clearly  seen  at  the 
single  frequency  of  the  sine  function.  In  the  application  to  the  gear 
mesh  induced  vibration.  Equation  (40)  indicates  that  n^  should  be  set 

equal  to  about  a  hundred  times  the  minimum  number  of  cycles  to  complete  a 
base  period  of  tooth  deflection  variation.  Therefore,  for  periodic  tooth 
deflection  variation,  Fourier  analysis  (subroutine  FOUR)  is  more  economical 
to  use  to  extract  the  frequency  content  of  the  induced  vibration,  since  it 
requires  the  tooth  deflection  data  over  only  the  minimum  number  of  mesh 
cycles  to  cover  one  base  period.  However,  for  random  variation  of  gear 
parameters,  only  subroutine  SPECT  can  be  used  to  obtain  the  frequency 
spectrum  of  the  induced  vibration. 
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APPENDIX  II 

COMPUTER  PROGRAM  FOR  PREDICTION  OF  GEAR  MESH  EXCITATION  SPECTRA 


Input  Variables.  Format,  and  Instructions 
Card  1  Title,  columns  2  through  72. 
Card  2  Control  numbers.  Format  (715) 


a.  NMC  Number  of  mesh  cycles. 

Place  the  last  digit  of  this  number  in  column  5. 

b.  INT  Identification  as  to  whether  this  control  card  rep 

resents  the  last  complete  set  of  input  data  being 
submitted . 

If  more  sets  of  input  data  follow,  use  0. 

If  this  is  the  last  set,  use  1. 


Plac2  this  digit  in  column  10. 


c.  MN  Classification  of  the  types  of  spur  gears  to  be  con¬ 
sidered  . 

If  both  the  driving  and  the  driven  gears  are  exter¬ 
nal  gears,  use  1. 

If  the  driving  gear  is  an  external  gear  and  if  the 
driven  gear  is  an  internal  (ring)  gear,  use  0. 

(The  program  will  not  run  properly  if  the  internal 
gear  is  submitted  as  the  driving  gear.) 

Place  this  digit  in  column  15. 


d.  MMM  Number  of  initial  terms  of  the  Fourier  analysis  for 
which  coefficients  will  be  printed,  beyond  the  co¬ 
efficient  for  the  constant  term.  This  number  cannot 
exceed  (FI  x  NMC  +  NMC/2),  where  FI  is  the  input 
variable  submitted  on  card  4. 


Place  the  last  digit  of  this  number  in  column  20. 

e.  IPLT  Instruction  as  to  whether  the  calculated  tooth 

meshing  error  is  to  be  plotted. 

If  tooth  meshing  error  is  to  be  plotted,  use  1. 

If  plotting  is  to  be  bypassed,  use  0. 

Place  this  digit  in  column  25. 

f.  IF0UR  Instruction  as  to  whether  Fourier  analysis  of  the 

tooth  meshing  error  is  to  be  performed. 

If  it  is  to  be  performed,  use  1. 

If  it  is  to  be  bypassed,  use  0. 


Place  this  digit  in  column  30. 


g.  ISPECT  Instruction  as  to  whether  power  spectral  density 
function  for  the  tooth  meshing  error  is  to  be 
calculated . 

If  it  is  to  be  calculated,  use  1. 

If  it  is  to  be  bypassed,  use  0. 

Card  3  Gear  design  data.  Format  (6E  13.5) 

a.  FN1  Number  of  teeth  in  the  driving  gear. 

Use  columns  1  through  13.  (Do  not  omit  decimal 
point .) 

b.  FN2  Number  of  teeth  in  the  driven  gear. 

Use  columns  14  through  26  .  (Do  not  omit  decimal 
point .) 

c.  RBI  Base  circle  radius  of  driving  gear,  in. 

Use  columns  27  through  39. 

d.  R01  Radius  to  the  outside  diameter  of  the  driving  gear, 

in.  This  should  be  reduced  by  any  radial  loss  in 
working  surface  at  the  tip  of  the  teeth,  as  from 
tip  rounding  or  chamfering. 

Use  columns  40  through  52. 

e.  R 02  Radius  to  the  outside  diameter  of  the  driven  gear, 

or  R12  if  external,  and  to  the  inside  diameter,  if  inter¬ 

nal,  in.  This  should  be  corrected  for  any  radial 
loss  in  working  surface  at  the  tip  of  the  teeth, 

as  from  tip  rounding  or  chamfering.  In  the  case  of 
an  internal  gear,  this  radius  must  be  equal  to  or 
greater  than  the  base  circle  radius.  No  check  for 
this  is  provided. 

Use  columns  53  through  65. 

Card  4  Gear  design  data,  continued.  Format  (6E  13.5) 

a.  RT1  Radius  to  the  beginning  (near  the  base  of  the  tooth) 

of  the  involute  profile  on  the  driving  gear,  in. 

This  is  used  in  the  program  only  in  a  design  check 
as  to  whether  adequate  length  of  involute  has  been 
provided  for  contact  on  the  teeth  of  the  mating 
gear  up  to  its  tip.  If  this  radius  is  not  speci¬ 
fied  in  the  gear  design  data,  this  check  may  be  by¬ 
passed  by  substituting  the  root  circle  radius. 


Use  columns  1  through  13. 


b.  RT2 


Radius  to  the  beginning  (near  the  base  of  the  tooth) 
of  the  involute  profile  on  the  driven  gear,  in. 

See  above  for  substitute  when  not  specified. 

Use  columns  14  through  26. 

c.  RM1  Radius  to  the  root  circle  of  the  driving  gear,  in. 

If  the  radius  submitted  is  smaller  than  the  computed 
base  circle  radius,  this  is  noted  in  the  output,  and 
the  input  value  of  root  radius  is  used  at  some 
points  in  the  program.  If  the  root  radius  is  suf¬ 
ficiently  smaller  than  the  base  circle  radius  so 
that  the  root  fillet  center  lies  inside  the  base 
circle,  the  tooth  outline  between  the  base  circle 
and  the  fillet  is  assumed  to  be  a  radial  line  by 
the  program. 

Use  columns  27  through  39. 

d.  RM2  Radius  to  the  root  circle  of  the  driven  gear,  in. 

For  the  case  of  an  external  gear,  the  same  com¬ 
ments  as  above  apply. 

Use  columns  40  through  52. 

e.  FI  Number  which  indirectly  establishes  the  number  of 

calculation  points.  The  number  of  these  points  will 
equal  one  plus  twice  the  value  of  FI.  The  calcula¬ 
tion  points  may  be  viewed  as  selected  contact  points 
on  the  true  involute  profile,  extended  where  neces¬ 
sary.  These  contact  points  with  the  mating  involute 
are  associated  with  specific  angular  positions  taken 
by  the  gear  as  it  is  rotated,  where  the  angular 
positions  correspond  to  uniform  subdivisions  of  the 
tooth  spacing  angle.  A  greater  number  of  these 
points  will  give  more  closely  spaced  point-by-point 
output  data.  A  greater  number  will  also  give  more 
accurate  calculations  of  tooth  deflections  and 
Fourier  coefficients.  A  value  of  FI  equal  to  12 
giving  25  calculation  points  has  been  found  to  be 
convenient . 

Use  columns  52  through  65.  (Do  not  omit  decimal 
point .) 


f .  T1 


Circular  tooth  thickness  at  the  pitch  circle  of  the 
driving  gear,  in.  The  radius  of  the  pitch  circle  is 
as  defined  in  card  3.  If  not  specified  in  the  gear 
design  data,  it  may  be  estimated  as  one-half  of  the 
difference  between  the  actual  circular  pitch  and  the 
working  backlash. 


Use  columns  66  through  78. 


54 


Card  5  Gear  design  data,  continued.  Format  (5E  13.5) 


a.  T2  Circular  tooth  thickness  at  the  pitch  circle  of  the 

driven  gear,  in.  The  comments  for  T1  also  apply 
here . 

Use  columns  1  through  13. 

b.  FI  Effective  tooth  face  width  of  the  driving  gear,  in. 

Where  the  face  widths  of  the  two  gears  are  similar, 
use  the  actual  face  width  without  any  reduction  for 
normal  end  chamfering  or  rounding.  Where  one  tooth 
is  much  wider,  use  as  its  effective  face  width  an 
amount  suitably  larger  than  the  narrower  width  to 
allow  for  the  limited  additional  support  that  the 
greater  width  provides. 

Use  columns  14  through  26. 

c.  F2  Effective  tooth  face  width  of  the  driven  gear,  in. 

The  comments  for  FI  also  apply  here. 

Use  columns  27  through  39. 

d.  RF1  Fillet  radius  on  the  driving  gear,  in. 

Use  columns  40  through  52. 

e.  RF2  Fillet  radius  on  the  driven  gear,  in. 

Use  columns  53  through  65. 

Card  6  Gear  material  properties.  Format  (6E  13.5) 

a.  YE1  Young's  modulus  (in  bending)  for  the  material  of 

the  driving  gear,  lb/in. ^ 

Use  columns  1  through  13. 

b.  YE2  Young's  modulus  (in  bending)  for  the  material  of 

the  driven  gear,  lb/in.^ 

Use  columns  14  through  26. 

c.  GE1  Shear  modulus  for  the  material  of  the  driving  gear, 

lb/ in. 

Use  columns  27  through  39. 

d.  GE2  Shear  modulus  for  the  material  of  the  driven  gear, 

lb/ in.^ 

Use  columns  40  through  52. 
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e.  P0S1  Poisson's  ratio  for  the  material  of  the  driving 

gear.  Since  this  ratio  is  used  only  in  the  allow¬ 
ance  for  the  "wide  beam  effect,"  it  should  be  re¬ 
duced  for  the  cases  where  tooth  face  width  is  not 
much  greater  than  tooth  thickness,  with  a  limiting 
value  of  zero  when  the  teeth  have  a  width  smaller 
than  the  thickness. 

Use  columns  53  through  65. 

f.  P0S2  Poisson's  ratio  for  the  material  of  the  driven 

gear.  Comments  for  P0S1  also  apply  here. 

Use  columns  66  through  78. 

There  are  NMC  sets  of  the  following  Card  7  and  Card  8.  Each  set  supplies 
center  distance,  tooth  spacing  errors,  tooth  profile  errors,  tooth  support 
compliances  and  tangential  load  for  one  mesh  cycle. 

Card  7  Center  distance  and  tooth  spacing  error  data.  Format  (3E  13.5) 

a.  CL  Center  distance,  in. 

This  must  be  the  actual  center  distance,  including 
any  substantial  spreading  under  load. 

Use  columns  1  through  13. 

b.  VPT1  Tooth  spacing  error  on  the  driving  gear,  in.  This 

error  is  based  on  the  distance  between  the  pitch 
points  of  successive  teeth,  but  the  error  is  ad¬ 
justed  to  apply  to  the  direction  of  the  line  of 
action.  This  adjustment  is  accomplished  by  multi¬ 
plying  the  pitch  line  error  by  the  cosine  of  the 
pressure  angle.  The  error  is  positive  if  the 
measured  spacing  is  smaller  than  the  desired 
spac ing . 

Use  columns  14  through  26. 

c.  VPT2  Tooth  spacing  error  on  the  driven  gear,  in.  The 

comments  under  VPT1  also  apply  here. 

Use  columns  27  through  39. 

Cards  8-1  to  8-2N  Point-by-point  data.  Format  (5E  13.5) 

Total  number  of  cards  equal  to  twice  the  number  of 
calculation  points  (Nj)  between  pitch  points  of 
adjacent  teeth,  or  the  same  as  two  plus  four  times 
the  value  of  FI  (see  card  4).  This  specifies  that 
cards  must  be  introduced  even  if  it  is  known  that 
there  is  no  contact  at  the  particular  calculation 
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a.  ZJ1 


point  or  even  if  the  tooth  profile  does  not  actual¬ 
ly  extend  to  the  calculation  point.  As  explained 
below,  a  blank  card  may  be  used  for  these  points. 

For  the  driving  gear,  the  first  card  is  for  the  cal¬ 
culation  point  located  (Nj-1)  points  preceding 
the  pitch  point  (or  inside  the  pitch  circle); 
the  (Nj)th  card  is  for  the  pitch  point;  the  last  or 
(2Nr)th  card  is  for  the  calculation  point  located 
(Nj)  points  after  the  pitch  point  (or  outside  the 
pitch  circle).  The  last  point  may  also  be  de¬ 
scribed  as  the  point  of  contact  on  one  meshing 
tooth  when  the  pitch  point  is  the  point  of  contact 
on  the  next  meshing  tooth. 

For  the  driven  gear  which  is  an  external  gear,  the 
first  card  is  for  the  calculation  point  located  (Nj) 
points  before  the  pitch  point  (or  inside  the  pitch 
circle);  the  (Nj+l)th  card  is  for  the  pitch  point; 
the  last  or  (2Nj)th  card  is  for  the  calculation 
point  located  (Nj-1)  points  after  the  pitch  point 
(or  outside  the  pitch  circle).  The  point  for  the 
first  card  may  also  be  described  as  the  point  of 
contact  on  one  meshing  tooth  when  the  pitch  point  is 
the  point  of  contact  on  the  previous  meshing  tooth. 

For  the  driven  gear  which  is  an  internal  gear,  the 
first  card  is  for  the  calculation  point  located  (Nj) 
points  following  the  pitch  point  (or  outside  the 
pitch  circle);  the  (Nj+l)th  card  is  for  the  pitch 
point;  the  last  or  (2Nj)th  card  is  for  the  calcula¬ 
tion  point  located  (Nj-1)  points  before  the  pitch 
point  (or  inside  the  pitch  circle).  The  point  for 
the  first  card  may  also  be  described  as  the  point 
of  contact  on  the  meshing  tooth  when  the  pitch  point 
is  the  point  of  contact  on  the  next  meshing  tooth. 

Deviation  of  the  point  on  the  actual  tooth  profile 
on  the  driving  gear  from  the  true  involute  (as  de¬ 
fined  by  the  gear  design  data),  in.  This  true  in¬ 
volute  is  positioned  relative  to  the  actual  profile 
so  that  its  deviation  at  the  pitch  point  is  zero. 
Where  the  deviation  represents  material  added  to 
the  true  involute,  it  is  positive;  where  it  repre¬ 
sents  material  subtracted,  it  is  negative.  The 
deviation  is  measured  normal  to  the  involute  pro¬ 
file.  If  the  profile  does  not  extend  to  the  partic¬ 
ular  calculation  point  or  if  it  is  known  that  the 
mating  gear  will  not  contact  at  this  point,  the 
deviation  may  be  noted  as  zero. 

Use  columns  1  through  13. 
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b.  UJ1 


Tooth  support  compliance,  or  any  compliance  supple¬ 
mentary  to  the  tooth  compliance  included  in  the 
analysis,  on  the  driving  gear,  in. /lb.  This  com¬ 
pliance  is  the  deflection  under  unit  load  at  the 
calculation  point  on  the  profile  in  the  direction 
of  the  load  (or  normal  to  the  profile).  A  uniform 
compliance  for  all  calculation  points,  such  as 
would  result  from  a  uniform  gear  shaft  compliance, 
would  not  affect  the  final  results  as  far  as  motion 
irregularities  or  load  transfer  is  concerned;  it 
would  only  increase  the  mean  deviation  in  transmit¬ 
ted  motion. 

Use  columns  14  through  26. 

c.  ZJ2  Deviation  of  the  point  on  the  actual  tooth  profile 

on  the  driven  gear  from  the  true  involute,  in.  The 
comments  under  ZJ1  also  apply  here. 

Use  columns  27  through  39. 

d.  UJ2  Tooth  support  compliance,  etc.,  on  the  driven  gear, 

in. /lb.  The  comments  under  UJ1  also  apply  here. 

Use  columns  40  through  52. 

e.  WT  Total  load,  tangent  to  the  pitch  circle,  transmitted 

by  the  gear  teeth,  lb. 

In  the  first  N  card  8s,  WT  should  be  left  blank. 

WT  in  the  second  N  card  8s  represents  the  loads  at 
the  N  calculation  points  in  one  gear  mesh. 

Use  columns  53  through  65. 

Card  9  Gear  Speed.  Format  (15, E  13.5) 

a.  NWS  Identification  as  to  whether  the  input  speed  is  the 

speed  of  driving  or  driven  gear. 

If  driving  gear  speed  is  inputted,  use  1. 

If  driven  gear  speed  is  inputted,  use  2. 

Place  this  digit  in  column  5. 

b.  WS  Driving  or  driven  gear  speed,  rpm. 

Use  columns  6  through  18. 

Output  Variables  and  Explanations 
1.  Title 


** 
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2.  Control  numbers  —  NMC,  INT,  NMZ,  MMM,  as  in  input  card  2. 

3.  Design  -  FI(-I),  FNl(-Nl),  RBI,  R01,  TR1,  RM1 

BLANK,  FN2(-N2),  BLANK,  R02  or  RI2 ,  RT2 ,  RM2 
Tl,  FI,  RF1,  YE1,  GE1,  P0S1 
T2 ,  F2 ,  RF2 ,  YE2 ,  GE2 ,  P0S2 
all  as  in  input  cards  3  through  6. 

Items  4-11  are  printed  for  every  mesh  cycle,  totally  NMC  cycles. 

4.  Mesh  cycle  identification  and  center  distance. 

5.  Input  listing  of  profile  error  and  supplementary  compliance  ~  ZJ1, 
UJ1,  ZJ2,  UJ2  as  in  input  card  8. 


6.  Pressure  angle,  degrees. 


7.  Incidental  data  -  RP1,  RBI,  BA1,  BR1,  ATI 

RP2 ,  RB2,  BA2,  BR2 ,  AT2 


where:  RPL 
RP2 
RBI 
RB2 
BA1 
BA2 

BR1 

BR2 


pitch  circle  radius  of  driving  gear,  in. 

pitch  circle  radius  of  driven  gear,  in. 

base  circle  radius  of  driving  gear,  in. 

base  circle  radius  of  driven  gear,  in. 

arc  of  approach  of  driving  gear,  rad. 

arc  of  approach  of  driven  gear,  rad  (negative  on 
internal  gears). 

arc  of  recess  of  driven  gear,  rad. 

arc  of  recess  of  driving  gear,  rad  (negative  on 
internal  gears). 


ATI  angle  of  rotation  of  driving  gear  from  the  position 
at  which  the  line  of  action  intersects  the  involute 
at  the  start  of  the  involute  profile  to  the  position 
at  which  the  line  of  action  intersects  the  involute 
at  the  pitch  point,  rad. 

AT2  similar  angle  of  rotation  of  driven  gear,  rad. 


Check  statement  when  part  of  the  profile  extends  within  the  base 
circle . 


Program  will  continue  in  any  case,  and,  where  necessary,  the  root 
radius  will  be  set  equal  to  the  base  circle  radius.  However,  in 
calculating  the  tooth  profile  and  the  tooth  deflections,  the 
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original  root  circle  radius  will  be  used  with  the  specified  fillet 
radius.  If  the  root  circle  lies  inside  the  base  circle  by  more 
than  this  fillet  radius,  a  radial  line  is  assumed  to  connect  fillet 
and  involute. 


8.  Driving  gear  data  -  Jl,  CJ1,  AJ1,  QJ1ABC,  XJ1,  YJ1,  XME1(=X), 
YMEl(-Y),  Jl,  QJ1A,  QJ1B,  QJ1C 


where:  Jl 

CJ1 

AJ1 

XJ1 

and 

YJ1 


XME1 

and 

YME1 

QJ1A 

QJ1B 

QJ1C 


identification  number  for  calculation  points  (see 
under  FI  of  card  4  in  the  input  data).  Listed  for 
values  of  (-21)  to  (21+1). 

condition  of  engagement  if  equal  to  one  and  no 
engagement  if  equal  to  zero. 

angle  of  rotation  from  the  position  of  contact  at 
the  pitch  point  to  the  position  of  contact  at  the 
calculation  point  -  negative  for  points  inside  the 
pitch  circle,  rad. 

coordinates  of  the  calculation  point  on  the  involute 
profile  with  the  origin  at  the  gear  center  and  with 
the  X-axis  as  the  centerline  of  the  tooth,  given 
only  for  the  points  at  which  contact  will  take  place 
with  the  mating  gear,  in. 

coordinates  of  the  point  on  the  root  circle  midway 
between  the  tangent  point  of  the  fillet  radius  and 
the  involute  profile  extended  (and  radial  inside 
the  base  circle),  in.  This  point  is  considered  to 
be  the  end  of  the  effective  base  of  the  tooth  for 
deflection  purposes. 

elastic  compliance  of  the  gear  tooth  acting  as  a 
cantilever  beam  in  bending  only,  normal  to  the 
profile  at  the  calculation  point,  in. /lb. 

elastic  compliance  of  the  gear  tooth  as  a  canti¬ 
lever  beam  in  shear  only;  otherwise  as  above. 

elastic  compliance  of  the  gear  tooth  as  a  rigid 
member  rotating  in  its  supporting  structure; 
otherwise  as  above. 


QJ1ABC  combined  compliance  of  the  three  above,  in. /lb. 

9.  Driven  gear  data  -  J2,  BLANK,  AJ2,  QJ2ABC,  XJ2,  YJ2,  XME2 ,  YME2 , 
J2,  QJ2A,  QJ2B,  QJ2C 

where:  J2  identification  number  for  the  calculation  points. 

For  external  gears,  J2  is  listed  for  values  of 
(-21-1)  to  (21).  For  this  case,  contact  takes 
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place  between  points  of  the  two  gears  for  which 
J1  =  -J2.  For  internal  gears,  J2  is  listed  for 
values  of  (21+1)  to  (-21).  For  this  case,  contact 
takes  place  between  points  of  the  two  gears  for 
which  J1  =  J2. 


All  other  variables  are  similar  to  their  counterparts  for  the 
driving  gear. 

10.  Input  tooth  spacing  error  data  —  VPT1,  VPT2 

as  in  input  card  7. 

11.  Tooth  meshing  errors,  loads  and  contact  compliance  —  JC1,  AJC1, 
EJT,  WTC ,  WTD,  WN,  WT,  QJD 


where:  JC1  identification  number  for  the  calculation  point  on 
the  first  tooth  of  the  driving  gear,  starting  with 
the  first  point  after  the  pitch  point  and  ending 
with  the  point  corresponding  to  the  pitch  point  of 
the  next  tooth. 


AJC1  angle  of  rotation  of  the  driving  gear  from  the 
position  with  contact  at  the  pitch  point  of  the 
first  mesh  cycle  to  the  position  with  contact  at 
the  calculation  point,  rad.  The  last  angle  in  the 
first  mesh  cycle  is  the  tooth  spacing  angle. 

EJT  tooth  meshing  error  or  deviation  from  pure  conju¬ 
gate  action,  as  a  pitch-line  linear  measurement  of 
the  motion  of  the  driven  gear  leading  the  driving 
gear,  in.  A  negative  value  indicates  that  the 
driven  gear  is  lagging  the  driving  gear,  as  might 
be  caused  by  deflection  of  the  teeth. 

WTC  tangential  load  carried  by  the  first  pair  of  teeth, 
lb. 


WTD  tangential  load  carried  by  the  second  pair  of 
teeth,  lb. 

WN  total  normal  load  transmitted  by  the  teeth,  lb. 

WT  input  tangential  tooth  load,  lb. 

QJD  contact  or  Hertzian  compliance  combined  for  both 
teeth  at  the  contact  point,  in. /lb. 


12.  List  of  tooth  meshing  error  over  1JMC  mesh  cycles  —  JC1,  AJC1,  EJT 

where:  JC1  identification  number  for  the  calculation  point. 

The  last  value  should  be  equal  to  (NMCxN) . 
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AJCl  same  as  AJC1  in  item  11. 
EJT  same  as  GJT  in  item  11. 


13.  Plot  of  tooth  meshing  error. 

Appears  only  if  IPLT  in  input  card  2  is  set  to  be  1. 

14.  Fourier  coefficients  —  I,  A,  B,  C,  KM 


where:  1 


A 


B 


C 


order  of  the  harmonic  to  which  the  coefficients 
apply.  The  zero  order  refers  to  the  constant 
component . 

the  Fourier  coefficient  of  the  consine  or  real 
component  for  that  harmonic  of  the  meshing  error, 
in.  The  value  for  I  =  0  is  twice  the  constant 
component  or  mean  value  of  the  meshing  error. 

the  Fourier  coefficient  of  the  sine  or  imaginary 
component  for  that  harmonic  of  the  meshing  error, 
in. 


o  o 

square  root  of  (A  +  B  ) .  Appears  only  if  IFOUR 
in  input  card  2  is  set  to  be  1. 


KM  ratio  of  frequency  to  tooth  mesh  frequency. 

The  following  output  is  concerned  with  the  power  spectral  density  function 
of  the  tooth  meshing  error.  Appears  only  if  ISPECT  in  input  card  2  is  set 
to  be  1. 


15.  Mesh  frequency,  cps 

16.  Incidental  data  “  FC,  FO,  BE,  H,  TR,  TMAX,  SM 

where:  FC  cutoff  frequency,  cps. 

FO  fundamental  frequency,  cps. 

BE  equivalent  resolution  bandwidth,  cps. 
H  sampling  interval,  cps. 

TR  total  record  time,  sec. 

TMAX  maximum  displacement,  sec. 

SM  maximum  lag  number. 

17.  Power  spectral  density  function  —  K,  FR,  GK,  GKK,  KM 

where:  K  harmonic  number,  from  0  to  SM. 


w 

v 
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frequency,  cps. 


'V 


FR 

GK 

GKK 

KM 


2 

"raw"  power  spectral  density  function.  In. 
"smooth"  power  spectral  density  function,  in.^ 
ratio  of  frequency  to  tooth  mesh  frequency. 
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'V 


PROGRAM  GGE AR (INPUT .OUTPUT . TAPES* I NPUT . TAPE6*0UTPUT .TAPE221 

•gear 

2 

014003 _ 

_ COMMON/ Y 7 /UJCL( 253  . U JD1 ( 25) . U J2 ( 25 )  .UJC2I25)  .UUl>*(25)  .ZJU50) 

Cll 6 

.  .J 

1  .ZJC1  <251 .2 J01  (25»  .2 J2 (501  iZJC2  (25)  .2 J02 (25)  .OJl  1501 ,0JC1 125)  .OJD1 

ggear 

4 

2(25)  .OJ2I50)  .0JC2I25) .  FJ1  (50)  .UJH50) .  CJC1(2  GGEAR 

5 

. 

— 3SI.CJOU25).  AJ1  (50)  «AJ2  (50)  .OJC  (50)  •  QJ0(  90)  .2X1501*230  tSO.  GGEAR 

6 

4)  «VJ(50)  »EJ1  (56)  .FJ2I50)  .AJC1  (25)  .AUC2(25)  ,FJC1  (27)  >0JO2(2S)  .ZM2( 

ggear 

7 

550) «UM2(50) .VTCI27) 

ggear 

8 

110803 _ 

_ COMMQN.PJ(501.CJ(5Q)  .OJ  (50)  tAJISO)  »AJ(50)  . YJI50) ) 

GGEAR 

_ 8. 

1A(15«0)  .1(1 500). 0(3000) 

MAY11 

1 

NHI1 

COMMON  D0.8R<BA.FJ.CC.YM«NN.N<FNJ«EP<TAN.RB.6M.F  tXM.MN.  II.M 

AL6 

1 

UI1I3 _ 

_  _  COMMON  (1A(50)  .06(50)  .OC  ISO) «QA1  (50)  t OBI  (50)  .OCI  (90>  <QA2(5fl).iQB2(S0  GGEAR 

_  12 

1)  .002(50)  .IH.IP 

GGEAR 

13 

•10003 

COMMON  NT  (SO). HEAD! (6) .HCA02I6) 

GGEAR 

14 

110003 _ 

_ COMMON/GU/NHCjINT »MMM. IPLT.FN1 . FN2.RB1  »R01.R02»Rf I.RT2.BM1  .RM2tFl.  AL0  . 

_ 2_ 

lTl.T2.Fl.F2.RFI.RF2.TEltYE2.GEl.GE2.P0Sl.P0S2.ANG 

ALG 

3 

•••M3 

DIMENSION  GGGI3000)  .WTCCI3000)  .FJCC1  (3000)  .AANGI3000) 

MAVll 

2 

uuu _ 

-  NR»5 

GGEAR 

. It. 

•00004 

NW«6 

GGEAR 

19 

•••005 

198  READINR.100)  (h£A01 1 1 ) • 1*1 .4) • (MEAD2 ( I ) • 1*1  .4) 

ggear 

20 

111121 _ 

_ -HEAD1MR.117L  NNC.INT.MN.MMK.IPLT.IFOUR.ISPECT 

MAY*. 

_ L. 

••••47 

REA0(NR,112IFN1.FN2.RB1.R01.R02 

GGEAR 

22 

••0065 

REAOINR.112IRT1.RT2.RM1.RM2.FI.T1 

ggear 

23 

M0105 _ 

READIMR.112)T2iFl  .F2.RF1.RF2 

ggear 

24. 

•00123 

READ (NR.  1 12)  YE1 .YE2.GE1 .GC2.P0S1.P0S2 

ggear 

25 

••0143 

WRITE (NW. 100)  (HE AD 1(1). 1*1.4). (HEAD? (I).I*1.4) 

ggear 

26 

•10163 _ 

_ WRITE.  INK.  1011 .  -  . 

GGEAR  . 

_ 21. 

•••167 

WRITE (NM.108I 

GGEAR 

28 

•••173 

WRITE (NW. 159)  NMC.INT.MN.MMM 

ggear 

29 

•••207 

WRITE  1NW.  104) 

ggear 

30 

••0213 

WRITE INW. 102)  F1.FN1.RB1.R01.RT1.RM1 

ggear 

31 

••0233 

IF IMN)  511.510.511 

GGEAR 

32 

M8234. - 

.  511-WRITE  INW. 1031  ...  . 

GGEAR. 

_ 33 

•••240 

GO  TO  512 

GGEAR 

34 

•••241 

511  WRITEINW.160) 

GGEAR 

35 

•••245 

512  WRITE  (MW.  158)  THETP.FN2.R02.RT2.RM2 

GGEAR 

36 

•00263 

WRITE (NW. 105) 

GGEAR 

37 

000267 

WRITE (NW. 102)  Tl. FI. RF1.YE1.GE1.POS1 

ggear 

38 

>00307 

- WRlXEiMM.107) . 

GGEAR  ... 

_ 30_ 

•00313 

WRITE (NW« 102)12 *F2»RF2*YE2»GE2.P0S2 

GGEAR 

48 

•00333 

WRITE (NW. 164) 

ggear 

41 

000337 

ANGwO.  . 

GGEAR 

42- 

••0340 

IMC-1. 

ggEar 

43 

••0342 

N*2.*FI«1. 

GGEAR 

44 

•10345 _ 

....  NN*N»N . 

ggear  . 

_ 45 

•••346 

1FIMN.EO.O)  T2*-T2 

JUNC14 

1 

••0351 

15*  WRITE (NW. 1500) 

MAY3 

2 

•10355 _ 

_  REAQ(NR.1121CL.¥PT1.VPT2 

ATTI 

1 

••0367 

WRITE (NW.1000)  INC. CL 

ggear 

47 

000377 

WRITE (NW. 1 1 0) 

final 

1 

000403 

...  t*NN _ 

CZIG 

2L 

•00405 

00  209  1*1 iNN 

ggear 

48 

•00406 

READ(NR.112)ZJ1(I).UJ1(I).2J2(I) . UJ2 (I).WT(I) 

GGEAR 

49 

•10423  ... 

. WRITE INW. 102) ZJl  ( I )  .UJ1  ( I ) .ZJ2 ( I ) .UJ2 ( I  > 

GGEAR 

50 

••0437 

2M2 (L>  *Z J2 ( I > 

GGEAR 

51 

64 
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RUN  VERSION  2.3  —  PSR  LEVEL  332—  GGEAR 


(00442 

UM2 (L 1 *U J2 ( I ) 

GGEAR 

52 

(00444 

IF(L-l)  200.209.225 

ggear 

53 

000446 

■  225  L*t-1 

GGEAR 

54 

0004SO 

209  CONTINUE 

GGEAR 

55 

0004S3 

CALL  6EAR0ICL.VPT1 .VPT2I 

GGEAR 

56 

000455 

NNN*INC*N 

GGEAR 

57 

(00460 

00  200  I«l.N 

ggear 

58 

000461 

JJ»NNN»I-N 

GU 

1 

0(0463 

AANGI JJt  *AJC1 (II 

MAY3 

3 

000466 

GGG< JJ)»EJI (I) 

GGEAR 

60 

0(0470 

G(JJ)*GGG(JJ) 

GGEAR 

61 

000472 

NTCC(JJ)*WTC(I) 

GGEAR 

62 

000474 

200  F  JCC1 <  JJt  * J J 

MAY2 

1 

000500 

IMC*IMC»1 

GGEAR 

64 

000501 

IF ( 1NC.LE .NHC )  GO  TO  ISO 

ggear 

65 

000503 

WRITE (NH. 1500) 

MAV3 

4 

0O0SO7 

WRITE  CNN. 2000) 

MAY3 

5 

000513 

DO  2500  1*1 .NNN 

MAY3 

6 

000515 

2500  WRITE (NW. 1021 F  JCC1  ( I )  .AANG ( I ) . GGG ( 1 1 

MAYS 

7 

000531 

IFIIPLT.EO.il  CALL  PLT  (FJCC1 .GGG.WTCC.NNNI 

GGEAR 

66 

000536 

M*F  I 

GU 

3 

000540 

IF  I IFOUR.NE.il  GO  TO  3000 

MAY9 

2 

000542 

WRITE (NW. 1281 

ALG 

4 

000546 

CALL  FOUR 

GGEAR 

67 

000S47 

LL*NNN»1 

ggear 

68 

(00551 

DO  250  1*1. LL 

ggear 

69 

000553 

LP*I-1 

ggear 

70 

000554 

SHK*1.*LP/NHC 

JUNE7 

1 

(00560 

RESLT*SORT ( A ( 1 1 ••2.BII  1  »#2I 

GGEAR 

71 

000566 

250  WRITE<NW.1?9ILP.A<I>  .81 1 1 .RESLT.SMK 

JUNE  7 

2 

000606 

3000  IFIISPECT.NE.il  GO  TO  3500 

MAY9 

3 

000610 

CALL  SPECT 

MAY9 

4 

000611 

3S00  CONTINUE 

MAY9 

5 

000611 

IFIINTI300.I98.300 

GGEAR 

73 

000612 

300  CALL  EXIT 

ggear 

74 

000613 

100  FORMAT  KA10I 

GGEAR 

75 

(00613 

101  FORMATI72HO  PN0335*COMPUTAT IONS  OF 

GEAR 

TOOTH  MESHING 

ERRORS 

GGEAR 

76 

I  3-22-1966  1 

ggear 

77 

000613 

102  FORMAT I7E13.51 

MAY9 

6 

000613 

103  FORMAT  (12H0  PRE.ANGLE7X2HN224X3HR1210X3HRT210X3HRM2I 

GGEAR 

79 

000613 

104  FORMAT  I/7X1MI 1 1X2HN1 1 1 X3HRB1 10X3HR01 1 0X3HRT  1 10X3MRM1 I 

ggear 

00 

000613 

1 05  FORMAT  1  /6X2HT 1 1 1 X2HF 1 1 1 X3HRF15X 1 2HY0UNGS 

MOD-I2X01HSHEAR  M0D-12X11 

ggear 

81 

1HPOS. RATIO-11 

ggear 

82 

000613 

107  FORMAT  (/6X2HT21 1  X2HF211X3HRF25X 12MY0UNGS 

M00-22XI 1HSHEAR  MOD-22XU 

GGEAR 

83 

1HPOS. RAT  10-21 

GGEAR 

84 

000613 

108  FORMAT I50H0  MESM  CYCLES  INPUT. 

DIFF 

.GEAR  NtOF  HARMONIC) 

GGEAR 

85 

000613 

109  F0RMATI4XI4. 4(8X141 1 

GGEAR 

86 

000613 

110  FORMAT I60H0INPUT  LISTING  OF  PROFILE 

ERROR 

ANO  SUPPLEMENTARY  COMPLI 

FINAL 

2 

1 ANCE/6X2H2 1 1 1 X2HU 1 1 1 X2HZ2 1 1 X2HU2 1 

final 

3 

000613 

112  FORMAT (6E13.5I 

ggear 

87 

0006L3 

117  FORMAT (8ISI 

MAY-9 

7 

000613 

128  FORMAT (43H0CALCULATED  FOURIER  COEFFICIENTS  FOR  E8RORS//7H  I9X  ALG 

5 

1SHAIII  13X4HBI II 14X1HC14X2HKMI 

JUNE7 

3 

(00613 

129  FORMAT (1 7. 3 C4X  E14.71.FU.7I 

JUNE? 

4 

000613 

158  FORMAT (2E13.5.13X.3E13.51 

GGEAR 

90 

010613 

160  -FORMAT ( 12H0  7X2HN224X3MR0210X3HRT210X3HRM21 

ggear 

91 

000613 

164  FORMAT I//16H  CALCULATED  OATAI 

GGEAR 

92 

(006 L3 

1000 .  F ORMATi  16H .  MESM .  CYCLE  NO..,  I3.20H. 

CENTER  OlSIONCEa, 

E13.61 . 

MAY  14.  _ 

JL 

0(0613 

1500  FORMAT (/6X  20(H*3XI/I 

MAY  3 

8 

000613 

2000  FORMAT (SX  3HJC19X4HAJC16X1 IHTAnG,  ERPORl 

MAY3 

? 

OOOM3 

_ STOP 

ggear 

94 

0(0615 

ENO 

GGEAR 

95 
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SUBROUTINE  GEARO(CL.VPTl.VPT2)  JUN13  1 

480046 _ CflMMON/Yl/iUCl  (2S)  iDJDl 125)  >UJ2  (25)  <U3C2(2S)  .  UJ02125)  .231150)  CZ1G  ..  .  3 


1.Z3C1 425)  .2  J01  (25)  .Z32 (50) .Z3C2(25>  <7302 (25)  .031(50)  ,QJC1  (25)  .0301  GEARO 

4 

2(25) .032(50) .03C2I25) •  F31 (50) .U3l(50) • 

C3CK2  GEARO 

5 

_ _  _ 

35) .C301 (25) .  A31  (50)  .A32 (50) »QJC (50) .030 ISO) .Z3C (50)  • 

2 3D (SO  GEARO 

6 

4) ,V3(50> .EJ1 (56) .F32(50)  .A3C1 (25) .A3C2I25) ,FJ Cl (27) .0302125) 

.ZM2I  GEARO 

7 

550) .UM2I50) »wTC (27) 

GEARO 

8 

-MO 044 - 

- COMMON  P3 (50) .03(50)  .03(50)  .43(50) «X3 (50) .Y3IS0) . 

GEARO 

9 

1 A (1500). B( 1500). G (3000) 

may!  i 

3 

000006 

COMMON  00<BR<BA<F 3.CC. YM.NN.N.FN3.EP.TAN.RB.GM.F .XM.MN. I I.M 

ALG 

7 

000006 

COMMON  04(50) .08(50)  .QC  (50)  .041 (50) .081 (50) .0C1 (50) .042(50) • 

082(50  GEARO 

12 

1).OC2(50)  . I M , I P 

GEAHO 

13 

000406 

COMMON  MT(50>.HE401 16)  .HEAD2I6) 

GEARO 

14 

000004  -  - 

C0MM0N/6U/NMC.  INT.MMM.  IPLT.FN]  .FN2.RB1  .ROl  .R02.HT  1  .RT2.RM)  .RM2.FI .  ALG 

B 

1T1.T2.F1.F2.RF1.RF2.YE1.YE2.GE1.GE2.POS1.POS2.ANO 

ALG 

9 

C 

SPECIFY  AND  INITIALIZE  READING  AND  WRITING  UNITS  FQR  IBM  1800 

GEARO 

17 

000006 

NR«5 

GEARO 

18 

004006 

NW»6 

GEAHO 

19 

000010 

198  CONTINUE 

GEARO 

20 

000014 - 

•  -  RRUCL«FN1/(FN1.FN2) 

APR30 

-  1 

000013 

IF (MN.EO.O) RP1«CL*FN1/ (FN2-FNI 1 

3UNE.7 

5 

000017 

CTT-RB1/RP1 

GEARO 

22 

M0021 

TTT«SQRT(1.-CTT»«2)/CTT 

GEARO 

23 

000026 

THETP-ATAMTTT) 

3UNE5 

1 

000030 

THETPP»THETP»1B0./3. 1415927 

3UNES 

2 

000033  - 

WRITE INM.112) 

YUNG 

1 

000036 

MRITE(NM.102I  TMETPP 

3UNE5 

3 

000044 

IFIF1-F2)  320.320.321 

GEARO 

25 

000451 

321  FF2*F2 

GEARO 

26 

000053 

GO  TO  322 

GEAHO 

27 

000053 

320  FF2»FI 

GEARO 

28 

4000S5 

322  -FF2*FE2**0.8 

GEARO 

29 

000061 

304  001— T1/2.0/RPLTHETP 

GEARO 

33 

000065 

F12»FN1/FN2 

GEARO 

34 

000067 

RP2*RP1/F 12 

GEARO 

35 

000070 

002— T2/2.0/RP2.THETP 

GEAHO 

36 

000074 

FN3«2.0*FI«1.0 

GEARO 

37 

-M0077 

N-FN3 

GEAHO 

38 

000101 

NN»N»N 

GEARO 

39 

000102 

CS-CTT 

3UNE6 

1 

000104 

TAN-TTT 

3UNE5 

4 

000105 

EP1»1.0-POS1»POS1 

GEARO 

43 

000107 

EP1-YE1/EP1 

GEARO 

44 

M0111 

EP2«I.O-POS2*POS2 

GEAHO 

45 

000114 

EP2«YE2/EP2 

GEAHO 

46 

000115 

EP12*( 1.0/EP1 *1  .0/EP2) *0.5 

GEARO 

47 

000121 

EP12-EP12»»0.9 

GEAHO 

48 

400125 

RB2-RB1/F12 

3UNE5 

5 

000127 

C1»R01/RB1 

GEARO 

51 

000130 

C2*R02/RB2 

GEARO 

52 

000132 

BR2» (SORT  (C2»C2-1.0)-TAN) 

GEARO 

53 

000140 

BR1« (SORT  (C1»C1-1.0)-TAN> 

GEAHO 

54 

000146 

C3*TAN/F12 

GEARO 

55 

000147 

IFIBR1-C3)  403.403.404 

GEAHO 

56 

N* 
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000154 

404  BR1*C3 

GEARO 

57 

000156 

WHITE  INK. 1261 BH1 

GEARO 

58 

000163 

403  BA1*8R2/F12 

GEARO 

59 

000165 

BA2*8R1*F12 

GEARO 

60 

000166 

IF(HN)  440.450.440 

GEARO 

61 

000172 

450  BA1—8A] 

GEARO 

62 

000173 

BA2*-BA2 

GEARO 

63 

000175 

440  IF  (BA1-TAM  401.401.402 

GEARO 

64 

000200 

402  BA1»TAN 

GEAtO 

65 

000202 

WRITE  INK. 125) BA1 

I  .ARU 

66 

000207 

401  C 1*6. 2831 845/FN1 

GEARO 

67 

000211 

IF(BAl-Cl)  213.213.611 

OtARO 

68 

000216 

611  WRITE  INK.  1 74) 

GEARO 

69 

000222 

213  IFIBRl'Cl)  215.612.612 

GEARO 

7>» 

000227 

612  WRITEINW.175) 

GEARO 

71 

000233 

215  C1*RT1/RB1 

GEARO 

72 

009235 

C2--RT2/RB2 

GEARO 

73 

000237 

Cl  -SORT  (C1»C1-1 .0)-TAN 

GEARO 

74 

000245 

AT1*ABS  (Cl) 

GEARO 

75 

000246 

C2  *50RT  (C2»C2-1.0) -TAN 

GEARO 

76 

000255 

AT2*ABS  <C2> 

GEARO 

77 

000257 

WRITE (NW. 122) 

GEARO 

78 

000262 

WRITE INW. 102) RP1. RBI. BA  1.RR1. ATI 

FINAL 

6 

000300 

WRITE  INK. 162) 

GEARO 

80 

000304 

WRITE! Nr  .102) RP2.RB2.BA2.0R2.AT2 

GEARO 

81 

000322 

NM/*) 

JUN13 

2 

000323 

NUZ*1 

JUN13 

3 

000324 

IF  CAT1-BA1 )  311.216.216 

GEARO 

84 

000331 

311  WRITE (Nr. 1 76) 

GEARO 

85 

000335 

216  IF IAT2-8A2)  309.217.217 

GEARO 

86 

000342 

309  IFIHN)  620.621.620 

GEARO 

87 

000343 

621  WRITE (NW. 1 78) 

GEARO 

88 

000347 

GO  TO  217 

GEARO 

8V 

000352 

620  WRITE (N«. 1 79) 

GEARO 

90 

000356 

217  IF(RMl-RBl)  191,192.192 

GEARO 

91 

000363 

191  AM1*1.0 

GEARO 

92 

000365 

WRITE (NW . 1 1 8 1 

GEARO 

93 

000370 

GO  TO  1920 

GEARO 

94 

000373 

192  AM1*RM1/RB1 

GEARO 

95 

000375 

1920  C1*AM| 

GEARO 

96 

000376 

AM1*SQRT  (AN1»ANI-1.0I-TAN 

GEARO 

97 

000404 

C2*RF1/R81 

GEARO 

98 

XT 

000406 

PF 1*C2 

GEARO 

99 

< 

000407 

C1«RM1/RB1 

GEARO 

100 

4 

000411 

C3= (Cl »C2) •*2-| .0 

GEARO 

101 

P 

'  000414 

IFIC3I  406.406.407 

GEARO 

i  02 

4 

000420 

406  C3*0.0 

GEARO 

103 

000421 

WRITEINW.156) 

GEARO 

104 

00042S 

TP01*0.0 

GEARO 

105 

! 

000426 

WRITE (NW. 156) POl 

GEARO 

106 

000434 

GO  TO  408 

GEARO 

107 

000437 

407  POl'SORT  (C3) 

GEARO 

108 

000441 

TP01*ATAN  (POl) 

GEARO 

109 

000444 

408  IF (RM2-RB2)  193. 194.194 

GEARO 

110 

000451 

193  AM2*1 .0 

GEARO 

111 

67 
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GEARO 


090*53 

WR!TE(NW,119> 

GEARO 

112 

000*56 

GO  TO  19*0 

GEARO 

113 

000*61 

19*_ AH2*RH2/RB2 

GEABQ 

...  11* 

000*63 

1960  Cl  MM2 

GEARO 

115 

000*6* 

AM2-SQRT  (AH2»AM2-1.0>-TAn 

GEARO 

116 

000*72 

C2-RF2/RB2 

GEARO 

117 

000*7* 

PF2* C2 

GEARO 

110 

000*75 

C 1 *RH2/RB2 

GEARO 

119 

-000*71.-..- 

-  1FIMN1  *56. *55. *56 

GEARO 

120.. 

000502 

*55  C2— C2 

GEARO 

121 

000503 

*56  C3* (Cl *C2I **2-1.0 

GEARO 

122 

900S06 

IF  (C3)  *09. *09, *)0 

GEARO 

123 

OOOSIO 

*09  C3*0.0 

GEARO 

12* 

000511 

WRITE (NW.1551 

GEARO 

125 

.400515 

TPD2*0.0 

GEARQ 

126 

000516 

WRITE <NW«1551PD2 

GEARO 

127 

00052* 

GO  TO  *11 

GEARO 

120 

000527 

*10  P02*S0RT  <C3) 

GEARO 

129 

000531 

TP02*ATAN  (P02) 

GEARO 

130 

OOOS3* 

*11  Kl»l 

GEARO 

131 

000535  .  - 

E1«RF1 

GEARO 

132 

000S36 

E2-PD1 

GEARO 

133 

0005*0 

E3«TP01 

GEARO 

13* 

0005*1 

E*»RM1 

GEARO 

13S 

0005*3 

00*001 

GEARO 

136 

0005** 

C1*AH1 

GEARO 

137 

0005*6 

201  C2*DD*C1 

GEARO 

130 

090550 

Cl'ATAN  (TAN*C1I 

GEARO 

139 

00055* 

C*«C1-C2 

GEARO 

1*0 

000556 

C3*E3-E2-D0»TAN 

GEARO 

1*1 

000562 

IF (KI-1 )  *22, *22, *25 

GEARO 

1*2 

000567 

*25  IF(RN)  *22, *21, *22 

GEARO 

1*3 

000570 

*21.  C3“-C3 

GEARO 

1** 

000571 

C*»-C* 

GEARO 

1*5 

000573 

*22  C1*C3*E1 

GEARO 

1*6 

000575 

C*«(C1»C*)*0.5 

GEARO 

1*7 

000600 

V»E**SIN  (C*l 

GEARO 

1*0 

000602 

X*E**COS  (C*l 

GEARO 

1*9 

000605 

IF(Kl-I)  202,202,203 

GEARO 

150 

000612 

202  XM1*X 

GEARO 

151 

000613 

VH1*V 

GEARO 

152 

000615 

K1*KI«1 

GEARO 

153 

000617 

00*002 

GEARO 

15* 

000620 

C1*AM2 

GEARO 

155 

000621 

E1-PF2 

GEARO 

156 

000623 

E2»PO 2 

GEARO 

157 

00062* 

E3*TP02 

GEARO 

150 

000626 

E**RN2 

GEARO 

159 

000630 

GO  TO  201 

GEARO 

160 

000630 

203  XH2*X 

GEARO 

161 

000631 

YN2*Y 

GEARO 

162 

000633 

EE*6.2831S5*/FNJ 

GEARO 

163 

000635 

cc*ee/fni 

GEARO 

16* 

000636 

FJ  *-FNJ»1.0 

GEARO 

165 

0006*0 

RB*RB1 

GEARO 

166 
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000641 

BR*BRI 

GEAHO 

167 

000643 

BAxBAl 

GEAMO 

168 

000644 

00*001 

GEAHO 

169 

000646 

XH*XMl 

GEAHO 

170 

000647 

YM*YM1 

GEAHO 

171 

000651 

GN*GE1 

GEAHO 

172 

000652 

F*F1 

GEARO 

173 

000654 

11*1 

CEARO 

174 

000655 

EP*EP1 

JEAHO 

175 

000657 

HRI TE  <NW* 151 1 

GEARO 

176 

000662 

CALL  AJCOH 

GEAHO 

177 

000663 

00  221  1*1 *NN 

GEARO 

178 

000667 

Cl  *XJ(I1 

GEAHO 

179 

000670 

C2  =  Y  J  ( 1 1 

GEAHO 

180 

000672 

F  J1 ( I ) *PJ III 

GEARO 

181 

000674 

C3*CJ( I ) 

GEAHO 

182 

000675 

AJH1)*AJ(11 

GEARO 

183 

000677 

CJ1(1)*0J(1) 

GEAHO 

1 B4 

000701 

IH1 a | H 

GEAHO 

185 

000702 

IP1*IP 

GEAHO 

186 

000704 

0A1 1 1 ) *0A ( I > 

GEAHO 

137 

000706 

OBI ( 1 1 *0R(  1 1 

GEAHO 

188 

000710 

0C1 II)*OC(l) 

GEAHO 

189 

000712 

IF (C3)  515.515,516 

GEAHO 

190 

000714 

515  MRITE(NH.10?)FJ1 (II .C3.AJ1 (1) 

GEARO 

191 

000726 

GO  TO  221 

GEAHO 

192 

000731 

516  WRITE(NW.102IFJ|  (II  ,C3.AJl(I).0Jt  (D.C1.C2 

GEAHO 

193 

000751 

221  CONTINUE 

GEAHO 

194 

000756 

HRITE(N6.120IXM1.tM1 

GEAHO 

195 

000765 

8RITE INN. 161  * 

GEAHO 

196 

000771 

00  517  I*IH.1P1 

GEAHO 

197 

000775 

517  MRITE(NH.102I*J1  (II  .OA1  (II.OBKII.OCl  (II 

GEAHO 

198 

001015 

00*002 

GEAHO 

199 

001016 

R8*RR? 

GEARO 

200 

001017 

xm*xm2 

GEAHO 

201 

001021 

YM*VM2 

GEAHO 

202 

001022 

IF(MNI  302.301,302 

GEAHO 

203 

001024 

301  FJ*FnJ 

GEAHO 

204 

001026 

GO  TO  310 

GEAHO 

205 

001026 

302  FJ*-FNJ 

GEAHO 

206 

001030 

310  F*F2 

GEARO 

207 

001031 

EP*E P2 

GEARO 

208 

001033 

C-H*GE2 

GEARO 

209 

001034 

CC*EE/FN2 

GEAHO 

210 

001036 

11*2 

GEAHO 

211 

001040 

VR I TE (NN , 1521 

GEAHO 

212 

001043 

CALL  AJCOH 

GEAHO 

213 

001044 

L*NN 

GEAHO 

216 

001046 

00  222  1*1 ,NN 

GEAHO 

215 

001051 

C4  *XJ(II 

GEAHO 

216 

001052 

C5  *Y J  ( 1 1 

GEARO 

217 

001054 

IH2*IH 

GEAHO 

218 

001055 

IP2*IP 

GEAHO 

219 

001057 

0*2(11 *QA (II 

GEAHO 

220 

001061 

QB2 ( 1 1 *OB ( 1 1 

GEAHO 

221 

69 
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GEARO 

•01063 

OC2U)*QCU) 

GEARO 

222 

•01065 

FJ2(!)*PJU> 

GEARO 

223 

_O01067 _ 

_ CltAJIIL 

fiFARO 

??4 

••1070 

AJ2ILUC1 

GEARO 

225 

••1073 

C2*0J< I ) 

GEARO 

226 

••1074 

0J2(L)«C2 

GEARO 

227 

•01076 

IF  (C2)  519.518.519 

GEARO 

228 

•01076 

518  WRITE(NW.132)FJ2UI.C1 

GEARO 

229 

OQUQA _ 

_ GO  10-520  ...  . 

..  HEARD _ 

730 

••1111 

519  WRITE (NW. 132) F  J2 (I ) .Cl .C2.C4. C5 

geaAo 

231 

••1127 

520  IFIL-1)  222.222.226 

GEARO 

232 

-Ml  133  _ 

226  L»L-1 

GEAAn 

233_ 

•01135 

222  CONTINUE 

GEAftO 

234 

•01140 

WRITE  INN. 120) XM2.TM2 

GEARO 

23S 

-101 147- 

_ WRITE  INM.  1631 

HEARD 

236. _ 

••1153 

00  521  I*IH2.1P2 

GEARO 

237 

••1157 

521  WRITE1NW. 102IFJ2I I) .042(1)  .082(1) »0C2(I> 

GEARO 

238 

-M1177 

00  251  1*1 .N 

GEARO 

239 

••1200 

IN*I»N 

GEARO 

240 

••1201 

FJC1  II)*FJ1  ( IN* 

GEARO 

241 

•01204 

- -AJCIUUAJIMN) 

GEARO.  _ 

.242 _ 

••1206 

AJC2U ) *A  J2 (IN) 

GEARO 

243 

••1210 

0J01(I)*0J1(I) 

GEARO 

244 

••1212 

0JC1  (t)*OJl  UN) 

GEARO 

245 

•01214 

0 JC2 ( I ) *0J2 (IN) 

GEARO 

246 

•01216 

0 JD2 ( 1 ) *0 J2 ( I ) 

GEARO 

247 

_ 

_  -XJOIUIpCJU) 

GEARO.. 

-248  _ 

••1222 

251  CJC1 II)*CJ  <IN) 

GEARO 

249 

••1226 

NUZ«I 

GEARO 

250 

M1227  - 

HNZ*1 

GEARO 

2S1 

••1230 

197  CONTINUE 

MAT  1 

1 

•01230 

L*NN 

GEARO 

253 

•01232 

-  -00-415  1*1. N 

GEARQ- 

254 

••1233 

IN*I«N 

GEARO 

255 

••1234 

UJ0KD-UJ1  U) 

GEARO 

256 

••1236 

UJC1  II)*UJ1  UN) 

GEARO 

257 

001241 

UJC2I I )  *UN2  UN) 

GEARO 

258 

001243 

U J02 ( I ) *UH2 U ) 

GEARO 

259 

•01243 - 

_ ZJ01U)*ZJ1U> 

GEARO 

260 

••1247 

ZJCI(I)-ZJIUN) 

GEARO 

261 

••1251 

ZJC2U )  *ZH2  UN) 

GEARO 

262 

••1253 

Z J02 ( I ) *ZH2 ( I ) 

GEARO 

263 

•01255 

0  JC  ( I )  *0 JC 1  U  )  «U JC 1  U  >  *0 JC2 1 1 )  *U JC2  <  I ) 

ge*Ro 

264 

001263 

0  JO  U )  *0  JO  1  U  )  *U  JO  1  II )  •  0  J02  ( I )  »U  J02 1 1 ) 

Gf  ARO 

265 

001271 

415  .ZJC  (I  >*ZJC1  ( I )  *ZJC2 1 1 ) 

GEARO 

266 

•01301 

19V  CONTINUE 

GEARO 

267 

••1301 

WR1TE1NW. 106) 

'»E  ARO 

268 

••1305 

WRITE INW. 102) VPT1.VPT2 

GEARO 

269 

•01317 

VP1*VPT1/CS 

GEARO 

270 

001322 

VP2*VPT2/CS 

GEARO 

271 

•01324 

418  CONTINUE 

MATl 

2 

•01324 

WRITEINW.115) 

GEARO 

275 

••1330 

00  501  1*1 .N 

GEARO 

276 

•01334 

WTU)*WT(N«I) 

GEARO 

277 

001337 

WN*WT(I)/CS 

GEARO 

278 

70 


GEARO 
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001341 

IFIHNI214.214.e00 

GEARO 

279 

001343 

800  CONTINUE 

GEARO 

280 

00134? 

CUHN*».l 

MAY1 

3 

00134? 

QD*1 ,37/FF  2*EP1 2/Cl 

MAY1 

4 

00I3S2 

ZJ' < 11*2 J01 < I )  »ZJD2 ( 1 1 «VP1-VP2 

GEARO 

281 

001360 

CCC  *0 JC ( 1 ) «QD 

GEARO 

282 

001362 

ODD  «OjD(II»QO 

GEAHO 

283 

001364 

C1*CJC1 1 1 1 

GEARO 

284 

001366 

C2«CJ01 (!) 

GEARO 

285 

001367 

IF (Cl  1  256«?56»257 

GEARO 

286 

001371 

256  IF (C2)  'IS. SOI. 258 

GEARO 

287 

001373 

257  IF (C?l  250.250.260 

GEARO 

288 

00I37S 

259  VJ(I!*?JC(II-CCC  *HN 

GEARO 

289 

001401 

HTC ( I i ~  f?N*C5 

GEARO 

290 

001403 

HTD*0.0 

GEARO 

291 

001404 

GO  TO  261 

GEARO 

292 

00140S 

258  VJ(II«ZJO(l)-OOD  *HN 

GEARO 

293 

001411 

HTC  ( 1 1 *0.0 

GEARO 

294 

001412 

HTO*HN*CS 

GEAHO 

295 

001414 

CO  TO  261 

GEARO 

296 

00141S 

260  C3*CCC  *HN‘2 JO < I ) -2 JC ( I ) 

GEARO 

297 

001422 

C4*0D0  *MN«ZJC(1I-ZJ0<I) 

GEARO 

298 

001426 

IF (C3I  259,259,262 

GEARO 

299 

001430 

262  IF (C4I  258.258,263 

GEARO 

300 

001432 

263  C1*CCC»000 

GEARO 

301 

001434 

C2*CCC  *230(11.000  *2JC(It-CCC*D0D*HN 

GEARO 

302 

001442 

VJIIXC2/C1 

GEARO 

303 

001444 

C1«CS/C1 

GEARO 

304 

00144S 

HTC( I1*ID00*HN.2JC ( I l-ZJO(I 1 1 *Cl 

GEARO 

305 

001452 

HTO* (CCC*HN.2 JO ( 1 1 -Z JC  ( 1 1 ) *C 1 

GEARO 

306 

001460 

261  EJ1(I)*VJ(I>/CS 

GEARO 

307 

001463 

AJC1  (D.ANG.AJCl  (I) 

GEARO 

308 

001465 

HRITE(NH.102IFJC1  (II  .AJCKII  ,EJ1  (1>«hTC(I).HT0.*N,iT(  n.OO 

MAY  1 

5 

001511 

501  CONTINUE 

GEARO 

310 

001516 

ANG«*JC1 (N) 

GEARO 

311 

001517 

24  IFINU2-MU2I  504,504,502 

GEARO 

312 

001522 

502  NUZ*NU2.1 

GEARO 

313 

001524 

CO  TO  199 

GEARO 

314 

001S24 

504  HUZ*I 

GEARO 

315 

001525 

IF (NM2-MM2 )  214,214.506 

GEARO 

316 

001530 

506  MM2*MM2.1 

GEARO 

317 

001532 

GO  TO  197 

GEARO 

318 

001S32 

214  CONTINUE 

GEARO 

319 

001532 

503  RETURN 

GEARO 

320 

001533 

102  FORMAT (8E13.5) 

MAY  1 

6 

001533 

106  FORMAT (//34H  INPUT  OATA  ON  TOOTH  SPACING  ERR0R/,6X4HVPT19X4MVPT2> 

GEAHO 

322 

001533 

110  FORMAT (60H0INPUT  LISTING  OF  PROFILE  FPROR  AND  SUPPLEMENTARY  COMPLI 
1ANCE/6X2H21 1 1X2HU1 1 1X2H221 1 X2MU2) 

GEARO 

GEARO 

323 

324 

001533 

112  FORMAT (16H0PRES. ANGLE  IOEG)  I 

FINAL 

7 

001533 

115  FORMAT  (42H0CALCULATEO  TOOTH  MESHING  tPPORS  AND  L0A0S//6X3HJC 1 9X4HA 

1 JC16X1  IhT ANG.  ERR0R6X3HHTCI 0X3HHTD 1 1 X2HHN1 IX2HMI 10X3HOJDI 

GEAHO 

MAY1 

325 

7 

001533 

118  FORMAT (74H0DRIVING  GEAR  INPUT  ROOT  RADIUS  SMALLER  THAN  BASE  CIRCLE 
1  RADIUS  1 

GEAHO 

GEARO 

327 

328 

00IS33 

119  FORMAT (74HOORIVEN  GEAR  INPUT  ROOT  RADIUS  SMALLER  THAN  BASE  CIRCLE 

1  RADIUS  1 

GEARO 

GEARO 

329 

330 
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GEARO 


001533 

.001533  _ 
001533 

001S33 

001533 
-001533  . 

001533 
001533 


001S33 

001533 

001533 

001533 

001533 

001S33 

001533 

00LS33 

001533 

001533 

001533 

001533 

001533 


120  FORMAT (48HOCOOO.  OF  EFFECTIVE  TOOTH  PROFILE  AT  MOOT  CIRCLEUX1HX13  GEARO 
1X1HY/52X.2E13.6)  GEARO 

_  122  FORMAT(/6X3HRP110X3HRB110X3HBAUOX3HBR110X3HAT1>  OLG  . 

125  FORMAT (66H0  DRIVEN  GEAR  TEETH  ENGAGE  UNDER  CUT  PORTION  OF  DRIVING  GEARO 

1GEAR  TEETH//6H  Bill".  E13.61  GEARO 

126  FORMAT I66H0  DRIVING  GEAR  TEETH  ENGAGE  UNDER  CUT  PORTION  OF  DRIVEN  GEARO 

1GEAR  TEETH//6H  BRlx  E13.61  GEARO 

132  FORMAT!  El 3.5 • 13X.4E 1 3.S )  GEARO 

151  FORMAT  I/6X2HJ1 1 1X3HCJ1 10X3HA  J18X6HQ  J1 ABC9X3HXJ1 10X3HY  Jll  GEARO 

152  FORMAT  (/6X2HJ224X3HA  J2SX6HOJ2ABC9X3HX  J210X3HYJ2)  GEARO 

155  FORMAT (63HOORIVEN  GEAR  INPUT  RADIUS  TO  FILLET  CENTER  INSIDE  BASE  GEARO 

1CIRCLE ./42HOPROGRAM  CONTINUES  HITH  CORRECT  TREATMENT. I  GEARO 

156  FORMAT (63H00RIVING  GEAR  INPUT  RAOIUS  TO  FILLET  CENTER  INSIDE  BASE  GEARO 

1CIRCLE./42NOPROGRAM  CONTINUES  WITH  CORRECT  TREATMENT.)  GEARO 

159  FORMAT (13X.4E13.S)  GEARO 

161  FORMAT  (/6X2HJ1 10X4HOJ1 A9X4HQJ1B9X4HQJ1C)  GEARO 

162  FORMAT (/6X3HRP210X3HRB210X3HBA210X3HBR210X3HAT2)  GEARO 

163  FORMAT  (/6X2MJ2 1 0X4HQ J2A9X4HQ J2B9X4H0  J2C )  GEARO 

170  FORMAT IN3M0CALCULATED  TOTAL  CONTACT  COMPLIANCE  0JD-.EI3.6)  GEARO 

171  FORMAT (23M0CALCULATEO  NORMAL  WNx  E13.6I  GEARO 

.  174  FORMAT  (SOHO  ANGLE  OF  APPROACH  ON  0R1VING  GEAR  IS  GREATER  THAN  TOOT  GEARO 

1H  SPACING  ANGLE.  PROGR AM/25HC0NT I NUE 0  WITHOUT  OVERLAP)  GEARO 

175  FORMAT  I75H0  ANGLE  OF  RECESS  ON  DRIVING  GEAR  IS  NQT  SMALLER  THAN  TO  GEARO 

10TH  SPACING  ANGLE.  /34H  PROGRAM  CONTINUED  MITHOUT  OVERLAP)  GEARO 

176  FORMAT  (80H0  DRIVING  GEAR  TEETH  MESHING  ON  PROFILE  INSIDE  OF  TIF  01  GEARO 

1AMETER.  PROGRAM  CONTINU-/2IHEO  WITHOUT  CORRECTION)  GEARO 


331 

332 
_  ID 

334 

335 

336 

337 
336 

339 

340 

341 

342 

343 

344 

345 

346 

347 
346 
349 
330 

351 

352 


355 

356 


178  FORMAT  (SOHO  ORIVEN  GEAR  TEETH  MESHING  ON  PROFILE  OUTSIDE  OF  TIF  D  GEARO  357 

1IAMETER.  PROGRAM  CONT IN-/22MUED  WITHOUT  CORRECTION)  GEARO  356 

179  FORMAT IBOMO  DRIVEN  GEAR  TEETH  MESHING  ON  PROFILE  INSIDE  OF  TIF  01  GEARO  359 

1AMETER.  PROGRAM  C0NTINU-/21HED  WITHOUT  CORRECTION)  GEARO  360 

END  GEARO  361 
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SUBROUTINE  AJCDH 

A  JCUH 

2 

•00802 

DIMENSION  SO  (SO )  .CD  (SO )  <8  IK  (SO  >  ,  A  IK  (5C  /  ,FLK  (SO)  .COCISO) 

AJCOH 

3 

•00002 

COMMON  PJ(50) .CJ<50) <QJ(50) < AJISO)  «XJ(50>  <YJ(50>  • 

MAY  1 1 

4 

1A ( ISO 0) <6 (ISOO) .0(3000) 

MAY  11 

S 

000002 

COMMON  DO.BR.BA.FJ.CC.YM.NN.N.FNJ.EP.TAN.RB.GM.F .XM.MN*  1 1  «’•• 

AJCUH 

s 

000002 

COMMON  OA  (SOI  .OB  (SO)  .QC(SO) .  0A1  (50)  .OBI  (50)  .OCl  iS0).QA2(5 

vB2 (50  AJCOH 

6 

1 ) .0C2 (SO  1 . IM. IP 

AJCOH 

7 

000  00 

00  106  I *  1 . NN 

AJCOH 

a 

••0CS4 

QA(I>*0.0 

AJCDH 

9 

000005 

06(11*0.0 

AJCOH 

10 

000006 

OF(I)“0.0 

AJCOH 

n 

000007 

0 J ( I >  >0. 0 

AJCOH 

12 

000010 

«J(l)*0.0 

AJCDH 

13 

000011 

106  YJ( I ) *0.0 

AJCOH 

14 

000016 

IF(II-l)  60S. 60S. 603 

AJCOH 

IS 

•00017 

60S  MM*1 .0 

AJCOH 

16 

•00021 

GO  TO  606 

AJCOH 

17 

000022 

403  NM*MN 

AJCDH 

18 

000026 

606  K*1 

AJCOH 

19 

•00025 

00  SOI  1*1. NN 

AJCOH 

20 

000027 

A J( I ) *F J*CC 

AJCOH 

21 

000031 

PJ( I ) *FJ 

AJCOH 

22 

•00032 

IF(lI-l)  612.612,413 

AJCOH 

23 

000035 

412  CALL  CALCJ(AJd)  .8R.0A .AVI 

AJCOH 

24 

000062 

CJ( I ) *XY 

AJCOH 

25 

000066 

GO  TO  414 

AJCOH 

26 

00006S 

413  XY*COC(I> 

AJCOH 

27 

000067 

414  IF (XV)  S03.S03.101 

AJCOH 

28 

000051 

101  IF(K-l)  102.102.103 

AJCOH 

29 

000056 

10?  IH*I 

AJCOH 

30 

000056 

1 01  X  X* A J ( I ) 

AJCOH 

31 

000060 

DJ*OD<AX 

AJCOH 

32 

000062 

SO ( I ) *S1N  (OJ) 

AJCOH 

33 

000066 

C0(I)*C0S  (OJ) 

AJCOH 

34 

000072 

C1*ATAN  (TAN. AX) 

AJCOH 

35 

000076 

C3*C0S  (Cl) 

AJCOH 

36 

000100 

C4*C1-0J 

AJCOH 

37 

000102 

IF (MM)  402,401.402 

AJCOH 

38 

000106 

401  C4*-C4 

AJCOH 

39 

000105 

402  C5*SIN  IC4I 

AJCOH 

40 

000107 

C6*C0S  (C4> 

AJCOH 

41 

0001 1 1 

C1*RB/C3 

AJCOH 

42 

000113 

X  J(  1 1 =C1 *C6 

AJCDH 

43 

000116 

V J( I 1 *C1 *C5 

AJCOH 

44 

000120 

IF(K-l)  242,242,243 

AJCOH 

45 

000122 

242  B1K(I)*(YJ(1) ••3» YM*»3) ?3.0*F 

AJCOH 

46 

000130 

AIK(l)*(YJ(II*YM)*F 

AJCOH 

47 

000133 

FLK(I)=XJ(I)-XM 

AJCOH 

48 

000136 

GO  TO  504 

AJCOH 

49 

000136 

243  BIMI>*(YJ(I)«»3.YJ(!-1>»3)/3.0*F 

AJC.UH 

50 

000166 

AIM|)*(YJ(I).VJ(I-1I)*F 

AJCOH 

51 

000167 

FLK(I)*XJ(I)-XJ(I-1) 

AJCOH 

52 

000152 

504  K*K.) 

AJCOH 

53 

000156 

503  IF(I-NN)  502.501,501 

AJCOH 

54 

73 


1  X 
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ajcoh 


•OOIST 

502  IF  (MM)  506.505.506 

AJCOH 

55 

•00160 

505  FJ*FJ-1. 

AJCOH 

56 

•00162 

60  TO  501 

ajCQh  . 

57 

000163 

506  FJ«FJ.1.0 

AjCOH 

SO 

000165 

501  CONTINUE 

AJCOH 

59 

•00170 

IF(II-l)  301.301.302 

AJCOH 

60 

•00172 

301  IF (MM)  601.302.601 

AJCOH 

61 

000173 

601  LL*NN 

AJCOH 

62 

•0017S 

00  303  1*1 .NN 

AJCOH 

61. 

•00176 

COC(I)»CJ(LLI 

AJCOH 

64 

000200 

IF (LL-1 )  303.303.305 

AjCOH 

65 

000203 

305  LL*LL-1 

AJCOH 

66 

000205 

303  CONTINUE 

AJCOH 

67 

•00210 

302  IP*K.lH-2 

AJCOH 

60 

000213 

410  00  246  L*IH.IP 

AJCOH 

69. 

000215 

£1*0.0 

AJCOH 

70 

•00216 

E2*0.0 

AJCOH 

71 

000216 

£3*0.0 

AJCOH 

72 

000217 

£4*0.0 

AjCOH 

73 

•00220 

EE*XJ(L)-XM 

AJCOH 

74 

•00223 

IM*L 

AJCOH 

15. 

000223 

IF (MM)  202.203.202 

AJCOH 

76 

000225 

203  EE*-EE 

AJCOH 

77 

000226 

202  DO  245  I=IH.IM 

AJCOH 

70 

•00230 

Cl  *XJ(LI-XJ(I> 

AJCOH 

79 

000233 

C5*FLM(I) 

AJCOH 

80 

•00236 

IF (MM)  207 .208.207 

AJCOH 

81 

•00236 

208  C1*-C1 

AjCOH 

82 

•00237 

C5«-C5 

AJCOH 

83 

•00241 

207  C2*C5/BIMI> 

AJCOH 

84 

000243 

C3*(C5.3.0*C1)*C5»3.0*C1»CI 

AJCOH 

85 

000251 

E4*E4«C2 

AJCOH 

86 

000253 

C4*2.0*C1«C5 

ajCdh 

87 

•00255 

E1*E1‘C3*C2 

AJCOH 

•8 

•00260 

E2»E2»C4*C2 

ajCdh 

•9 

000263 

245  E3*E3«C5/AIK(II 

AJCOH 

90 

000271 

IF (MM)  206.205.206 

AJCDH 

91 

000272 

205  IM*IM-1 

AJCOH 

92 

000274 

206  C1*C0(L) *CD (L) 

AjCOH 

93 

000276 

C2*SO(L)*YJ(L) 

AJCOH 

94 

000300 

OAILI  *C1/3.0*E1/EP. (C2*E4-C0(L) •E2)/EP*C2 

AJCOH 

95 

000312 

OB(L)  *C1/6M*E3*1  .2 

AJCOH 

96 

000317 

C1*(EE*C0<L)-YJ(L)»SD<L) )/YM 

AJCOH 

97 

000323 

0C(L)*C1/EP*1.327/F*CI 

ajCdh 

98 

000330 

246  QJ(L)*0A(L) *QB(L> «QC(L) 

AJCOH 

99 

000336 

RETURN 

AJCDH 

100 

000337 

END 

AjCOH 

101 

RUN  VERSTON  2.3  —  PSR  LEVEL  332 


"V 


SUBROUTINE  FOUR 

FOUR 

2 

C  REQUIRES  2*N» 1  POINTS  Fill 

FOUR 

3 

C  POINTS  CORRESPOND  TO  ThET  A*2*PI / (2*N. 1 ) < . . , »2*PI 

FOUR 

4 

C  OUTPUT  A  ( 1 ) «B ( I )  REFER  TO  COSINE  AND  SINE  OF  <1-1 » »THETA 

FOUR 

5 

000002 

COMMON  Pj  (50  •  »C  J(50  >  .CIJ  <50 )  «  AJ  (50 1  •  A  J  (50  »  .Y  J(50>  « 

FOUR 

6 

1 A  1 15001 >R 11500) .6(30001 

MATH 

6 

000002 

COMMON  OO.BR.BA.FJ.CC.YM.NN.N.FNJ.EP.TAN.RB.GM.F.XM.MN.II.M 

FOUR 

8 

000002 

COMMON/OU/NMC *  INT  .MMM.  1PLT  .FN1  .Fn2  «RB  l  »R01  .R02.RI  1  *RT2*RM1  .RM2  .F 1 • 

MAY9 

8 

1T1.T2.F1.F2.WF1.RF2.YE1.YE2.GE1.GE2.POSI.POS2.ANO 

MAY9 

9 

000002 

C=1.0 

FOUR 

9 

000003 

S*0.0 

FOUR 

10 

000004 

M  =  IFI*(Fr.«NMC»NMC/2 

MAY9 

10 

000010 

M«M.| 

FOUR 

11 

000012 

NA*N*NMC 

MAY9 

11 

000014 

N2*NA-1 

MAY9 

12 

000015 

T 1 *2./NA 

MAY9 

13 

000020 

T  2=11*3.1415427 

FOUR 

15 

000021 

ClrCOS  (T2) 

FOUR 

16 

000023 

S1*SIN  (T2I 

FOUR 

17 

000026 

00  7  IPal.M 

FOUR 

18 

000027 

U1*0.0 

FOUR 

19 

000030 

U2*0.0 

FOUR 

20 

000031 

DO  3  1=1. N2 

FOUR 

21 

000032 

J=N2-I«2 

FOUR 

22 

000034 

U3=G(J).2.0»C«U1-U2 

FOUR 

23 

000042 

U2=U1 

FOUR 

24 

000043 

3  U1*U3 

FOUR 

25 

000047 

A(IP)*T1*(G(1) »C*U1 «U2I 

FOUR 

26 

000054 

B(IP|aTl»S«Ul 

FOUR 

27 

000056 

AA=A(IP)«C-8(IP)«S 

FOUR 

28 

000062 

BB*A(IP)»S«B(1P»«C 

FOUR 

29 

000065 

A ( IPI =AA 

FOUR 

30 

000067 

0 ( IP) =BB 

FOUR 

31 

000071 

Q=C1*C*S1*S 

FOUR 

32 

000074 

S=C1»S.S1«C 

FOUR 

33 

000077 

7  C*0 

FOUR 

34 

000103 

RETURN 

FOUR 

35 

000103 

END 

FOUR 

36 

-SUBROUTINE  CALCJ  (Cl  .C2.C3.C4) 

CALCJL.  - 

_ 2  _ 

••••C? 

IF (C3)  240.231 >241 

CALCJ 

3 

NT 

••0011 

DATA  CHEA03/4H  LB./ 

CALCJ 

4 

V 

••0011 

241  IF (Cl I  239.218. 218 

CALCJ 

5 

A 

••0013 

230  IF (ABS  (CD-C3)  220.220.219 

CALCJ 

6 

f 

•00017 

231  IF (Cl  1  219.218.218 

calCj 

7 

4 

••0021 

2.40-1F1C11  219.219.242 

CALCJ  . 

. 

••••23 

242  C3»-C3 

calCj 

9 

••0024 

IFIC1-C3)  219.218.218 

calCj 

10 

; 

•00026 

216  1FCC2-C1)  219.220.220 

CALCJ 

11 

000030 

220  C4.1.0 

CALCJ 

12 

•00031 

GO  TO  221 

CALCJ 

13 

••0032 

-219  C4»0.0 

CALCJ 

14 

••••33 

221  RETURN 

calCj 

15 

••0034 

END 

CALCJ 

16 

75 


f  * 


( 


RUN  VERSION 

2. 

3  --PSR  LEVEL  332— 

SUBROUTINE  PLT(X,Y,Z.NPTS| 

PLT 

2 

000007 

DIMENSION  T ( 3000). X (3000) ,2(3000)  .XXX  14) .XYY (4) ,022 (4) • 1LABX (5) , 

NAY11 

7 

i  ILA8Y (5) .HEX03 (3) 

PLT 

4 

000007 

COMMON  PJ(50)»CJ(50).QJ<501 .  AJ(50) .X J(50)  .YJ(50> * 

MAY  11 

0 

14(1500). BU500). 0(3000) 

MAYll 

9 

000007 

COMMON  00,8R.8A,FJ,CC»YM,NN.N.FNJ,EP.TAN.R8»GM»F*XN,MN. H.M 

PLT 

6 

000007 

COMMON  04(50)  .08  (50 )  .QC  (50 )  .041  (50)  .081  (50).OC1  (50)  ,042(50)  .082(50 

PLT 

7 

1 ) ,002(50)  . IH, IP 

plt 

8 

000007 

COMMON  WT(SO) .MEA01 (6) .HE 402 (6) 

maTI 

11 

000007 

DATA  (ILABX(J)  ,J*1.3)/10H  C4LCUL,  IOH4TION  P0IN.2HT  / 

PLT 

10 

000007 

OATA  ( ILABY ( J> , J* 1 ,3) /I OMT ANGENT I4L , 1  OH  ERROR  (IN.2H.)/ 

plt 

11 

000007 

DATA  ILAB2/IOHLOAO  (LB.)/ 

PLT 

12 

000007 

OATA  <HEA03(l).l*l,2>/10HTANG.  FORC.amE  »  / 

PLT 

13 

000007 

IN*S 

PLt 

14 

000007 

10*6 

plt 

IS 

000010 

ICAL*22 

PLt 

16 

000012 

CALL  PLOTSIIBUF.4OO0.ICAL) 

PLT 

17 

000014 

CALL  PL0T(4.,4.5,-3) 

PLT 

IB 

000017 

CALL  PLOT (-2., -4. ,3) 

PLT 

19 

000022 

CALL  0ASHPT(-2.,6.,.2SI 

plT 

20 

000025 

CALL  SCALE (X.5. «NPTS»« 1 ) 

PLt 

21 

000032 

Y (NPTS*1 1*0. 

PLT 

22 

000036 

00  10  K*1.NPTS 

PLt 

23 

000040 

10 

Y(K)*-Y(K) 

plt 

24 

000043 

CALL  SCALE (Y«4.,NPTS»1«*1) 

plt 

25 

000050 

Y(NPTS.2)»Y(NPTS*3) 

PLf 

26 

000056 

CALL  SCALE (Z,3.,NPTS»»1) 

PLt 

27 

000061 

AXX 1 1 )*X (NPTSM ) 

PLT 

2B 

000065 

AXX(2I*X(NPTS»2) 

PLT 

29 

000067 

AYY(1)*V (NPTS* 1 ) 

PLT 

30 

000071 

AYY(2)*Y(NPTS*2) 

PLt 

31 

00007] 

A22())*2(NPTS.l) 

PLT 

32 

000074 

AZ2(2)*7(NPTS.2) 

PLT 

33 

000077 

CALL  AXIS(0.,0.,ILABX,-22.5.,0.,AXX(I )  ,AXX(2>) 

PLT 

34 

000106 

CALL  AXIS (0..0., ILABY. 22.4,. 90..0..AYY  (2) ) 

PLT 

35 

000116 

CALL  FLINE(X.Y.-NPTS. 1,1,3) 

PLT 

36 

000125 

CALL  PLOT(0..-4.,-3) 

PLt 

37 

000130 

CALL  AXIS (0. ,0. , ILABX, -22,  5. ,0. .AXX ( 1 1  ,AXX (2) ) 

PLt 

38 

000140 

CALL  AX!S(O..O.i ILA82, 1 0,3. ,90, , A2Z  ( 1 ) , 422(2) ) 

PLt 

39 

0001S0 

CALL  FLINE (X.2.-NPTS. 1,1,11) 

PLT 

40 

000157 

CALL  PLOT  < .5,9. 5, -3) 

PLT 

41 

000162 

CALL  SYMBOL  (0. ,0. »•  l4»HEADl .0. .40) 

PLT 

42 

000166 

CALL  SYMBOL  (0,  ,-.25,  •  14  .ME 402,0,  ,40) 

PLT 

43 

000172 

CALL  SYMBOL (0.,-.5».14, HE 403,0., 14) 

PLt 

44 

000176 

CALL  NUMBER (999. ,999... 14, *T.O., 2) 

PLT 

45 

000202 

CALL  SYMBOL (999., 999.,. 14.EHEAD3.0., 4) 

PLT 

46 

000206 

CALL  PLOT (6. ,-9. ,-3) 

PLT 

47 

000211 

CALL  OAShPT (0.,9.,.25) 

PLt 

48 

000214 

CALL  PLOT (0. .0 • ,999) 

PLT 

49 

000217 

CO  20  K*1 .NPTS 

PLt 

50 

000223 

20 

Y (X) *-Y (K) 

PLT 

51 

000226 

RETURN 

plt 

52 

000227 

ENO 

plt 

53 

76 


i 


run  version 


£00002 

000002 

000002 

000002 

000002 

000003 


00000* 

0000*4 

400016 

000020 

000024 

000032 

000035 

000037 

000041 

000044 

000046 

000051 

000052 

000055 

000057 

000062 

000074 

000100 

000114 

000115 

000117 

000123 

000125 

000127 

000132 

000134 

000136 

000137 

000140 

000144 

000147 

000151 

000153 

000157 

000165 

000172 

000174 

000175 

000177 

000206 

000207 

000211 


2.3  —  PSR  LEVEL  332— 

»tr^r“^^“,,E2*°Ei,Ge*,P0S1 

Wl*5 

»:!&!?  5  BSMBW- 

RE  AD (NR • 100)  NVS«VS 
WS*WS*6. 28/60. 

FN.FNl*NS/6.28  28 

IF  (NWS. EO. 2)  fh«nf 

WRITE (NWtllO)  ER 
FC>FM«N/2. 

NT*N»NMC 

TR«NMC/FM 

H.l./FH/N 

FO-l./TR 

M*NT/20 

SM*M 

TMAX*N*H 

eE*l./TMAx 

SKiffiS! 

SUM*0  « 

00  10  1*1  *F*T 
10  SUM*SUH»G(l) 

SA*SUM/N7 
00  20  1*1 »67 
CN(1)*G(I)”S* 

20  CONTINUE 
MM*M«1 

DO  30  I*1«RR 
1R«I-1 

FR(D*IR*EC/SM 

CC*1./*nT",r* 

RR ( I ) *0  * 

00  25  J*1*NT 

C  RR(1>*R(I"1*  nt>  60  To  30 
25  SR{|)*GN(j)*GN(J.IR»»«R‘>' 

30  RR(I)*RR<1'  CC 

ms*h-i 

00  40  K*1»Rh 

KK(K)*A-1  DOiylM,,_l. )•***(«) 

GK(K)*RR(1),bs,mr’  ’  l* 

SSS*°. 

DO  35  JK*1»ms 
JJ*JK*1 


MAT9 

hath 

hath 

MATH 

GEARC 

HATH 

MATll 

MATH 

MATH 

MATll 

MATll 

MATH 

MATll 

MATll 

math 

MATll 

MATll 

MATll 

MATll 

MATll 

MATll 

HATH 

MATH 

MATll 

MATH 

MATH 

MATll 

MATll 

MATll 

MATll 

MATll 

Kmt  1 1 

MATll 

MATH 

math 

MATll 

MATll 

MATH 

MATll 

MATll 

MATH 

MATll 

MATll 

MATll 

MATll 

MATll 

MATll 

MATll 

MATll 

MATll 

MATH 

MATll 

MATH 

MATH 


14 

10 

11 

12 

1 

H 

15 

16 

17 

18 

19 

20 
21 
2i 

23 

24 

25 

26 

27 

28 
?9 
33 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
2 
3 


77 


RUN  VERSION  2.3  — PSR  LEVEL  332—  SPEC! 


000212 

*A*3.1416*JK*KK(K)/SN 

MATH 

62 

000217 

35  SSS*SSS*2.»RR(JJI»COSIAA) 

MAY  1 1 

63 

000230 

40  GK(K)*(GK(K)*SSSI*2.#H 

MAY  1 1 

64 

000236 

GKK  <  1 1 *.5*GK<1 ) «.5*GK (2> 

MAY11 

65 

000241 

GKK(MMJ*.5*GK(M).  »5*GK  ( MMI 

MAY  1 1 

66 

000246 

DO  50  I.L-2.R 

MAY11 

67 

000247 

GKK  (LH  *  .25»GK  ILL-1 )  ♦  .S*GK  (LL »  ♦  ,?5«GK  (LL  ♦  1 1 

MAY  1 1 

66 

0002SS 

50  CONTINUE 

MAY1 1 

69 

000260 

WRITEINw.130! 

MATH 

70 

000264 

WRITECNW.I40I 

MAY  1 1 

71 

000270 

DO  60  1*1 .PR 

MAY  1 1 

72 

000272 

SKM*KK  (I)*N/(2.*SMI 

JUNE  7 

6 

000277 

WRITE(NW.1501KK(1).FR(I)  .GKtD.GKKU)  .SKM 

JUNE  7 

7 

600314 

60  CONTINUE 

MAY  1 1 

74 

000317 

100  FORHAT(I5.E13.5) 

MAY  1 1 

75 

000317 

110  FORMATI/1X26MMESM1NG  FREQUENCY  IN  CPS  * »E 1 3.51 

MAY  1 1 

76 

000317 

115  FORMAT  (/7A2HFC1 1 A2HF011X2HBE) 

MAY  1 1 

77 

000317 

120  FORMAT (BE  13. 5) 

MATH 

78 

000317 

125  FORMAT  »/7*lHHia*?MTRll*4MTMAX9*2MSM) 

MAYll 

79 

000317 

130  FORMAT  </l*31HP0*ER  SPECTRAL  DENSITY  FUNCTION! 

MAY  1 1 

80 

000317 

140  FORMAT </6*.lMK6*2HFRll*2MGKll*3HGKK10A2MKM) 

JUNE  7 

8 

000317 

150  FORMAT  (I7.3EI3.5tF13.5l 

JUNE  7 

9 

000317 

RETURN 

MAYll 

83 

000317 

END 

MAYll 

84 

78 


